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1 Introduction: what is an effective theory? 

An effective theory is a systematic approximation to some underlying dynamics 
(which may be known or unknown) that is valid in some specified regime. An 
effective theory is not a model, since its systematic character means that, in prin- 
ciple, predictions of arbitrary accuracy may be made. However, if this is to be true 
then a small parameter, such as the a of quantum electrodynamics, must govern 
the systematic approximation scheme. As we shall see here, in many modern effec- 
tive theories the expansion parameter is a ratio of two physical scales. For instance, 
in effective theories of supersymmetric physics "beyond the standard model" the 
ratio would be of p or m, the momentum or mass of a standard model particle, to 
MgusY- In an effective theory for finite-proton-size effects in the hydrogen atom the 
small parameter would be Tp, the proton size, divided by rb, the Bohr radius. The 
smallness of this parameter is then indicative of the domain of validity of the effec- 
tive theory (ET). In this sense effective theories, like revolutions, carry the seeds of 
their own destruction, since the failure of the expansion to converge is a signal to 
the user that he or she is pushing the theory beyond its limits. Within the radius of 
convergence of the ET the ET "works" because of the following fundamental tenet: 

Phenomena at low energies (or long wavelength) cannot probe de- 
tails of the high-energy (or short-distance) structure of particles. 

I suspect that most physicists subscribe to this tenet: indeed, if the tenet were not 
true, then physics (other than calculations using a "theory of everything" ) would 
be impossible. 

In this first section I will begin by discussing the basic ideas of effective theo- 
ries using a few simple examples from undergraduate physics. In this way we will 
move from classical effective theories, to classical effective field theories (EFTs), to 
quantum effective field theories. Of course, to understand the latter one must be 
able to compute quantum-field-theoretic loop graphs, so this requires a little more 
education than the standard undergraduate curriculum contains — at least the un- 
dergraduate curriculum at American Universities! Nevertheless, I will attempt to 



Czechoslovak Journal of Physics, Vol. 51 (2001), No. 



A 



1 



Daniel Phillips 



present the work at a level that is understandable by someone who has completed a 
first course on quantum field theory, but has not necessarily yet studied the theory 
of regularization and rcnormalization. 

In Section 2 I will begin to focus on the nucleon-nuclcon (NN) system. This will 
first necessitate some definitions of terms, notation, and so forth. I then move on 
to discuss the special issues stemming from the presence of shallow bound states in 
the NN problem. After displaying one solution to this difficulty, I will define and 
employ an effective field theory which takes into account the presence of shallow 
bound states, but other than this only contains neutrons and protons as dynamical 
degrees of freedom. I will give examples of the success of this EFT, known as 
EFT(7f), in computing (very-)low-energy electromagnetic observables in the NN 
system. 

The extension of this work to the three-body problem raises some intriguing 
problems of rcnormalization. I will attempt to elucidate these in Section 3, where 
I draw on the work of Bedaque, Hammer, and van Kolck, to show how, after some 
thought and interesting discoveries, EFT(7f) can be applied to the NNN problem. 

Finally, in Section 4 I give a brief tour of results in an effective field theory 
with pions. The effective field theory of QCD in which nucleons and pions are the 
degrees of freedom is chiral perturbation theory (xPT). We were fortunate to have 
one of the founders of xPT, and indeed a pioneer in the field of EFTs, lecturing at 
this school. Prof. Leutwyler's lectures in this volume should be read in conjunction 
with the material presented here. Indeed, this article should be regarded as little 
more than light reading on the subject of nuclear EFTs. It is very far from being a 
thorough review on the topic. The reader who wishes to study the subject in detail 
should consult the reviews Ref. ||] which contain much more information than 
does this manuscript. These two reviews also contain full references to the original 
literature, a job I have not tackled in any systematic way here. 



1.1 Some very simple effective theories 

1.1.1 Gravity for /i < i? 

One of the simplest effective theories I know is one that is learned by high-school 
physics students. It concerns the standard formula for the gravitational potential- 
energy difference of an object of mass m which is raised a height h above the Earth's 
surface: 



AU = mgh, (1) 

where g is the acceleration due to gravity 

Of course from Newton's Law of Universal Gravitation (itself an effective theory, 
valid in the limit of small space-time curvature), we have 

GMm GMm , , 

AC/- + , 2 

rf n 
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for an object whose distance from the Earth's centre is initially ri and finally rf. 
Now, if we write 

Tf = n + h, (3) 
and assume that rj « R, the radius of the Earth, then 

GM R , 

AC/ = m — h. (4) 

R^ R + h ^ ' 

Identifying = g we see that Eq. (Q) is only the first term in a series, which 
converges as long as h/R <1: 

AU = mgh(l-^ + ^ + ..)j . (5) 

If we try to apply this theory to a satellite in geosynchronous orbit (h 3> R) the 
series will not converge. But for the space shuttle {h ^ a few hundred km) this 
series should converge fairly rapidly. Equation (|l|) is the first term in the effective 
theory expansion ^ for the gravitational potential energy, with that ET being 
valid for h < R. 

1.1.2 Effective theories in the hydrogen atom 

Presumably, the hydrogen atom is ultimately described in terms of string theory, 
or some other fundamental theory of physics. Nevertheless, to very good precision, 
we can use quantum electrodynamics (QED) as an effective theory to compute 
its spectrum. The reason why we can ignore corrections to QED from physics at 
the Planck scale when calculating the hydrogen-atom spectrum will become clear 
below. 

In the case of the hydrogen atom there is an effective theory for QED that is 
valid up to corrections suppressed by one power of a = e^/(47r), the fine-structure 
constant. That effective theory is known as the Schrodinger equation with the 
Coulomb potential. The radial wave function Uni{r) obeys the differential equation 
0:5 

1 d'^Uni + a 

Unlir) Unl[r) = Er,Unl[r). (6) 



2mo dr^ 

The solution, for the lowest-energy eigenstate (n = 1; ^ = 0) is, of course: 

wio(f) — A/" exp(— amor) = N exp(— r/rb) , (7) 

where rb — (amc)"^ ~ (4 keV)~^ ~ 0.5 A, and M is determined by the normaliza- 
tion condition. The corresponding eigenvalue is 

1 fi'-' 



^ ) Throughout I work in units where h = c = 1. 
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The Bohr radius, rb, sets the scale for most phenomena associated with the electron 
in the Hydrogen atom. In particular, the typical momentum of the electron is 
~ 1/rb, which means that relativistic corrections to the Schrodinger equation are 
suppressed by (mcrb)^^ = a^, thereby validating the non-relativistic treatment. 
Note that if we were discussing muonic hydrogen the energy and distance scales 
would be very different, since « 2M''b- 

The Bohr radius is large compared to the size of the proton, rp, and also com- 
pared to the scale of internal structure of the electron. One sense in which the 
electron has internal structure is that it is dressed by virtual photons. In fact, the 
Lamb shift is, in fact, just such an electron-structure effect, so "finite-electron-size" 
effects must be considered if very accurate results are desired. If we for the moment 
ignore the electron's internal structure and consider only the internal structure of 
the proton we would replace the Coulomb potential —a/r, by the potential gener- 
ated by an extended proton: 

with p(r') the local electric charge density of the proton at the point r'. Now we 
make a multipole expansion, in order to obtain: 



^« = -^E(7) y dVp(r')P„(r.f')(y for 



(10) 



with the nth Legendre polynomial. Here, p(r') only has support for r' < rp, and 
so the integrals are all numbers of order one. Since the solution of the differential 
equation (^) is mainly sensitive to r ^ rb, the expansion parameter here is rp/rb ~ 
1 A/1 fm, and so this series converges rapidly, with it entirely permissible to evaluate 
the corrections for the finite size of the proton in perturbation theory. Nevertheless, 
an accurate computation requires inclusion of the term of order (rp/rb)^ in this 
expansion. 

In fact, this term of 0((rp/rb)^) is the first correction due to finite-size effects 
in Eq. (^). This is easily seen from Eq. (^0|), since the coefficient of the term of 
0(rp/rb) is zero, as long as the proton's charge distribution is even under parity. 
Thus consideration of the scales in the problem alone would lead us to grossly 
overestimate the magnitude of the finite-size effect. It is the combination of scales 
and symmetry that leads to an accurate estimate of the magnitude of the effects 
neglected by assuming that the proton is point-like in Eq. (^). These two principles: 

— a ratio of scales generating an expansion parameter, 

— symmetries constraining the types of corrections that can appear, 
inform the construction of most effective theories. 
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1.2 Building a (classical) effective theory: the scattering of light from atoms 

Next I want to discuss an example of ET-construction which first appeared in print 
in the effective field theory lecture notes of Kaplan Q . These lecture notes are an 
excellent introduction to EFT, and this example provides a great demonstration of 
ET construction. Here 1 have reworked some of the notation, but the basic idea is 
as in Ref. 

Consider the scattering of low-energy light from an atom. The question we must 
answer is: What is the Hamiltonian that describes the interaction of the atom with 
the electromagnetic field of the incoming light? To do this, we first consider the 
scales in the problem; the energy of the electromagnetic field, uj is assumed small 
compared to the spacing of the atomic levels, AE, and the inverse size of the atom. 
We will assume in turn that all of these scales are much smaller than the mass of 
the atom. Using rb to estimate AE we have the following hierarchy of scales: 

LU<S:AE ^ < — < A/atom. (H) 

n^c^b Th 

Meanwhile the symmetries of the theory will be electromagnetic gauge invari- 
ance [C/(l)om], rotational invariance, and Hermitian conjugation/Time reversal. 
These symmetries constrain the types of terms that we can write in our Hamilto- 
nian. Firstly, the constraint of gauge invariance suggests that it is wise to construct 
-ffatom out of the quantities E and B, rather than using the four- vector potential 
A^. Then, secondly, rotational invariance suggests that we employ quantities such 
as and E • B in i/atom- However, E • B is odd under time reversal, and so we 
cannot write down a term proportional to it in if atom- Meanwhile, terms such as 
V • B, and V x E may be included in i?atom, but then can be eliminated from the 
Hamiltonian using the field equations for the electromagnetic field in the region 
around the atom: 

V-E = 0; V-B = 0; (f2) 
VxE = 0; VxB = 0. (f3) 

This leaves us with: 

i/atom = aiB2 + asE^ + asidoB^) + a^idoE'') + a^iE • B)^ + . . . (14) 

In spite of our cleverness in constraining the terms that may appear, we are still 
left with infinitely many operators that can contribute to ffatom- How are we to 
organize all of these contributions? 

Interlude: naive dimensional analysis 

The answer lies in the scale hierarchy established in Eq. (pi]), together with a 
simple technique known as naive dimensional analysis (NDA). This works as follows: 
consider, for instance, the operator B^. Counting powers of energy/momentum we 
see that it carries four powers of energy, which we write as: 

[B']=4. (15) 
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However, 7?atom must have dimensions of energy. It follows that ai and 02 must 
each carry three negative powers of energy/momentum: 

[ai] = H = -3, (16) 

that is to say: 

ai, 02 = 7 r^- (17) 

(some energy scale) 

Now we ask what energy scales are present in the problem and so might appear 
in the denominator here. The photon energy lu cannot appear in the denominator 
since the scale that occurs there should refer to a property of the atom. Any of 
the scales r^^, AE, or A/atom could be involved though. The most conservative 
estimate would be that: 

"i'"^^(A^- 

However, very low-energy photons cannot probe the quantum level structure of 
the atom: they should interact with the entire atom in an essentially classical way. 
Thus, AE cannot occur in the denominator of this lowest-dimensional term in the 
Hamiltonian, and so we deduce that ai and 02 must scale with ro, i. e.: 

01,02 -^Tq, (19) 

where the ^ usually indicates that the coefficient here could be a 3 or a 1/3 (or a 
-3 or a -1/3) but will generally be a number of order one |^. ai and 02 are in fact 
proportional to the electric and magnetic polarizabilities of the atom Q| . 

Meanwhile, the operators multiplying the coefficients 03 and 04 have dimension 
5. Thus, [03] = [04] = —4, and so 03 and 04 carry one more energy-scale downstairs 
as compared to ai and 02. Conservatively, we assign the scaling: 

03,04 (20) 

Similar estimates can be made for the other terms in -ffatom- The key point is that 
since and are the lowest dimension operators allowed by the symmetries 
and not already constrained by field equations, they will give the dominant effect 
in -ffatom for low-energy photons. Any higher-order effects will be suppressed by at 
least oj/AE, i. e.: 



aiB^ + a2E2 + O (^) 



(21) 



where oi and 02 are now dimensionless numbers. It is straightforward to turn this 
result into a prediction for the photon-atom cross section. Since [a] = —2 and the 



^ ) There is a subtlety here: since this is an electromagnetic interaction the argument here suffices 
to get the scaling with ut correct, but it does not count powers of «cm = 1/137 which is an 
additional small parameter in the problem. 
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cross section results from squaring the quantum-mechanical amplitude arising from 
the Hamiltonian (|2l|) we discover that 



,4^6 



(22) 



where the ~ disguises the hard work needed to figure out all the factors of 2, tt 
and so forth that really go into deriving cr! The strong dependence of ct on a; is, of 
course, the reason the sky is blue — as pointed out in Ref. H or in Ref. d, where 
the constant of proportionality in Eq. (^) is worked out in detail! 

1.3 A classical effective field theory: Fermi electroweak theory 

Equation ( pl| ) is an effective expression for the classical Hamiltonian that is valid at 
long wavelength, or equivalently, for low-energy electromagnetic fields. In general 
effective field theories are derived for low energies, although this need not be the 
case. 

A canonical example of a low-energy effective field theory is Fermi's electroweak 
theory. This is an effective field theory that can be used to compute, say, low-energy 
electron-neutrino scattering. The only particles explicitly appearing in this theory 
are electrons and neutrinos. By contrast, in the standard model, the electrons and 
neutrinos interact by the exchange of W and Z bosons. If we wish to compute the 
scattering of neutrinos from electrons we could compute the full standard model 
amplitude for diagrams such as Fig. |l| §: 



with: 



A = — (vPLiuiy) o — (n'^Qe) ; (23) 

^^2cos6iwy g2_^2 \^2cos6'w/ 



Pl,r = i(lT75); (24) 

Q = (-l + 2sin2 6'w)^'L + 2sin2 6'w^'R. (25) 

and q = p' — p being the change in momentum of the neutrino. Relating to 
laboratory quantities, we see that: 



q^ = ^AE,^^EU,in'^^j, (26) 

where Ei^h {E^^^) and ^lab are the initial (final) energy and scattering angle of the 
neutrino in the lab. system. So, if i?iab <C A/z, then we can expand the Z-propagator 
in Taylor series. The leading term in this series is then: 

A = i\/2GF (i^PlTm^^) (e7^Qe) , (27) 

with: 

G = ^ ^ g" ^28) 

^ 4V2cos2 6'w ^'^1 4%/2M2,' 
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Fig. 1. A Feynman diagram for neutrino-electron scattering via exchange of a Z particle. 



The effective Lagrangian that generates the amplitude (|2j) is: 

^Fcrmi = Gf\/2 (z^7^PlJ^) {e^'Qe) ■ (29) 

This is the amplitude written down by Fermi many years ago for beta-decay and 
associated processes. The argument given here for deriving G from the underlying 
theory goes under the rubric "tree-level matching" . 

Note that while the range of the neutrino-electron interaction (|2^) is of order 
1/A/z, the interaction ( ^9|) is of zero range. All four fields in the Lagrangian (^ ) 
are to be evaluated at the same space-time point, x. This approximation is sensible 
provided that the neutrino's energy is such that its Compton wavelength is much 
longer than l/Mz- In this circumstance the standard model Lagrangian may be 
approximated by a Lagrangian in which the electrons and neutrinos interact only 
when they are right on top of each other, since then the neutrino cannot resolve 
the finite range of its interaction with the electron. A more accurate approximation 
to the standard-model may be obtained by expanding the Z-boson propagator in 
powers of <7^/M|. These higher-order pieces of the Fermi electroweak theory La- 
grangian will give rise to operators of a higher dimension than the operator ([29|), 
which has dimension: 

2[i/] + 2[e] = 6. (30) 

This is a crucial point, since the operator Dvee is the lowest-dimension one we 
can construct out of the four lepton fields. This means that even if we knew nothing 
about electroweak interactions, we would still expect an interaction like (^9|) to be 
the leading term describing low-energy electron-neutrino scattering. Furthermore, 
since [C] = 4, and [ip] — 3/2, we have [Gp] = —2, and thus if we also had some 
inkling that electroweak physics was determined by the one energy scale Mz ~ 
Mw ~ 100 GeV, we would expect: 

^ (100 Gcvy ^^^^ 

Of course, this argument is ahistorical, since Fermi did now know Mz or Afw- 
However, the argument can be turned around: if Fermi had suspected that beta 
decay was due to the exchange of some heavy particle he could have inferred its 
approximate mass from measurements of Gp. 
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1.4 Effective quantum field tiieories: Scalar theory in 5+1 dimensions 

An essential — perhaps the essential — feature of quantum field theories is that one 
can compute loop corrections to the tree- level result. Now in fact Fermi Electroweak 
Theory was used in the previous section purely as a classical field theory of Dirac 
particles, since only the tree-level amplitude (|2^) was discussed. We cannot be 
said to have derived an effective quantum field theory until we have explained 
how to compute loops. Here we encounter an apparent difficulty, because the ef- 
fective field theories we have written down in the previous two subsections are 
"non-renormalizable" quantum field theories. Interactions such as (|29|) are highly 
singular at short distances, and this leads to divergences when loop graphs are 
computed. Such divergences are not an insurmountable difficulty: they occur in 
quantum electrodynamics, which is rendered predictive by the use of subtractions 
in the misbehaving integrals and the input of two pieces of physical data (mo and 
e). However, in "non-renormalizable" field theories such as Fermi electroweak the- 
ory the calculation of higher-loop graphs leads to divergences with more and more 
complicated structures. So much so that at each new order in the loop expansion 
additional pieces of physical data are needed in order to render the theory finite. 
This appears to be a bad situation, since it suggests that the theory has no pre- 
dictive power, as it requires infinitely many inputs before it can make a prediction. 
This conclusion is, however, hasty. In this section I will explain how to compute 
loop graphs in a non-renormalizable effective field theory, and indicate why the 
theory has predictive power in spite of its non-renormalizability. 

This will be done using an example effective field theory which has no real uses 
that I am aware of. Although perhaps this has changed with the recent interest in 
dilatons, radions, and so on, in higher-dimensional field theories. The field theory 
is that of a massive but light scalar particle in 5+1 dimensions. We will take the 
symmetries of the theory to be six-dimensional Lorentz covariance and invariance 
under the discrete transformation cf) —^ —(f>. This leaves us with infinitely many 
interaction terms we can write down in the Lagrange density: 

c = ^(a,</-)(a^</.) - im^^^ - _ ^2 A ^4 _ _ C (32) 

Here [C] —6, [cf)] —2, and so [A] — —2, implying that all the interaction terms here 
are non-renormalizable. However, the statement [A] = —2 also implies that 

A=^, (33) 

where a is a number of order one, and we have assumed that this field theory is an 
EFT for an underlying theory with one key heavy scale, A. Likewise: 

r o, ^ b ^ c ,„,, 
'=A^-' ^=A^= ^-A^' (^4) 

where we will assume that a, 6, and c are all of order one, i. e. they are "natural" . 
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Fig. 2. Feynman diagrams for (fxf) scattering at threshold in six-dimensional (f)"^ theory. 
The first hne represents the leading-order result, from the vertex X(f>'^, while the second 
line includes the one-loop graph constructed from these leading-order vertices, as well as 
the next-to- leading order vertex m^A<^*. 



What would we then expect for the threshold amphtude for (p(j) scattering on 
the basis of NDA? In general, Tth will be a function of to, the light scalar mass, 
and so, since [Tth] = —2, we expect Tth to vary in the following way: 

1 TTl^ / 2 \ 

Tth = TTo + ^T77 + • ■ • + possible terms proportional to log ( -?W I , (35) 

here M is just a convenient mass scale determined by the size of the threshold 
amplitude at to = 0, and H is order one. 

By contrast, direct calculation with (^) shows that the "leading-order" (LO) 
contribution to Tth, coming from the graph on the first line of Fig. |^ is just: 

T^t'h^ = (36) 

But leading-order in what? In fact, it will prove to be leading order in an expansion 
in to/ A. And we already see that if the scale M is of order A, which it should be 
for a generic quantum mechanics problem, then adjusting a to reproduce the first 



term in Eq. (35) will only require that a be of order one. Systems where M is not 
of order A — such as the NN system — are referred to as "fine-tuned" and require 
special treatment. We will discuss this in the next section. 

The next order in the to/ A expansion is provided by the two graphs on the 
second line of Fig. |^. Both of these graphs are 0{m'^ /A"^). However, the loop graph 
has a logarithmic divergence multiplying to^/A^. Using dimensional regularization 
to evaluate the loop integral in 6 — e dimensions we find that the overall contribution 
of the two graphs to Tth is: 

= - Y (4^ ^ (I - ^) + 

(If you are unfamiliar with dimensional regularization you should consult a field 
theory book such as Ref. Q for a primer on this topic.) Making a Laurent expansion 
of the Gamma function for small e, and expanding the other functions present in 
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powers of e too, we find that as e — > 0: 



rpNLO 



1 



2 (47r)3 



log 



+ A 



(38) 



where c is an m-independent constant, whose computation is left as an exercise: 

Exercise 1 Show that c = 1 — 7 + log(47r), where 7 is the Euler-Mascharoni con- 
stant. 

Putting in our NDA expectations A = 02 /A^ and A = a/A^ we see that T^'^ 
is indeed of order w? j and so will be suppressed relative to T]^ as long as the 
scalar is light relative to A. However, it appears that T^'^ is a small infinity, since 
part of the expression for it diverges as e ^ 0! In order to resolve this dilemma we 
need a piece of physical data in order to fix the counterterm A. In just the same 
way, in QED the bare electron mass is chosen to be infinite, with that infinity used 
to cancel off infinite contributions from the one-loop self-energy, and a physical 
parameter, namely mc, must be specified in order to render the theory predictive. 

In our case we will suppose that we know the coefficient H in Eq. (|35|), from 
calculations in the underlying theory, say: 



th 



T' 



(39) 



In this case A must be chosen so that: 

A2 1 



r^ = A + 



2 (47r)3 



2 

- +c- 
e 



log ( 



(40) 



Note that the constant A is renormalization-scheme dependent: choices of "minimal 
subtraction" (MS) or "modified minimal subtraction" (MS), will result in different 
definitions for A. However, the physical prediction for Tth is independent of renor- 
malization scheme: if we invert ( |40| ) and use it to remove A from ( |38| ) then we find 
that the infinity as e — > has been absorbed in A, and: 



rpNLO 



A2 



1 



2 (47r)- 



1 [ ™o 
■ log — J 



(41) 



What we see here is that the coefficient of the logarithmic term that describes 
the non-analytic variation of T^'^ with m is 'predictive, and is calculable order-by- 
order in the low-energy theory. This non-analytic behaviour arises from the long- 
distance, low-momentum, infra-red physics of the light-particle loops. The short- 
distance physics is encoded in the coefficient A, or, equivalently, in the piece of 
physical data Tg. 

The example we have examined here is really quite simple, but the behaviour 
it displays is generic. If we consider any loop graph contributing to the process 
— > (^(\) then its superficial degree of divergence is: 



d — riidi 4- 2L, 



(42) 
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where Ui is the number of vertices of type i, di is the number of powers of momen- 
tum that this type of vertex carries and L is the number of loops in the graph. 
However, if dimensional regularization and a mass-independent subtraction scheme 
are employed the sum of this loop and appropriate counterterms must carry d 
powers of m and/or momentum p. Hence it must behave like: 



where / is a dimensionless function of p and m. This has the further consequence 
that in order to ensure that all of the appropriate counterterms are present, and so 
all divergences can be removed, we must include all possible local counterterms in C 
which scale as (p, m)*^. In other words, we must include in our original Lagrangian 
all possible terms of that order which are consistent with the symmetries of the 
theory. Ultimately this amounts to calculating the loop graph in order to obtain 
the function / and then making as many subtractions as are necessary to render 
/ finite. Thus the effective field theory, although non-renormalizable, can make 
predictions for the mass and energy dependence of observables. 

A key point here is that in any non-renormalizable theory the degree of diver- 
gence d of graphs increases as more vertices and more loops are added. This had 
been regarded as a "bug" in non-renormalizable field theories, but now we see it is 
instead a "feature" |^. Non-rcnormalizability becomes a "feature" because it means 
that graphs with more loops and/or more complicated vertices carry more powers 
of the small scales (p, TOtt); and so are suppressed relative to tree-level and low- loop- 
order graphs. This fascinating paradigm shift is discussed further in Refs. H, ^, p^ . 

In a non-renormalizable field theory computing all loop graphs and tree-level 
contributions that carry up to d powers of p and m,r we will obtain a prediction for 
the (f)(j) — > amplitude which is accurate up to terms which are suppressed by an 
amount p^ j Is? and m? / relative to those that are included. EFT predictions such 
as the one obtained above for the mass dependence of the scalar-scalar scattering 
amplitude in six dimensions can thus be systematically improved by computing 
more loops, including more counterterms etc. 

1.5 Summary: the EFT algorithm 

Thus we can develop what I shall refer to as "The EFT algorithm". The steps 
in this procedure, which can be used to build a low-energy theory of just about 
anything, are (with apologies to the experts who will recognize that this is somewhat 
simplified) as follows: 

1. Identify the degrees of freedom of interest. 

2. Identify the low-energy scales and high-energy scales in the theory. The ratios 
will form expansion parameters, and so a good scale separation will facilitate 
construction of the EFT. 

^) Credit to Tom Cohen for this remark. 
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3. Identify the symmetries of the low-energy theory. 

4. Choose the accuracy required. This, together with the size of the largest 
expansion parameter, say x, will determine the order n to which observables 
must be computed. 

5. Write down all possible local operators, consistent with the symmetries, which 
have dimensions up to that order: 



6. Derive the behaviour of loops, that is to say, the analog of Eq. ( [42| ) for this 
particular theory. This will tell you how many loops you need to calculate to 
have a prediction that is accurate to nth order in x. 

7. Calculate the loops and renormalize them. 

At this point you should have a prediction which will be accurate up to terms 
whose relative size will be set by x (as long as the mass scales which generated x 
can actually produce a term at the next order). Now, in fact, to check that the 
theory is really converging in the expected manner an n + 1st order calculation is 
necessary. Comparison with experimental data can also be made, but note that in 
addition to the usual experimental errors EFTs carry an intrinsic theoretical error, 
of size x'^'^^. This error is an honest reflection of the impact of our ignorance of 
high-energy physics on the observable in question. 



In this section I will implement this "EFT algorithm" for nucleon-nucleon scat- 
tering at low energies, i. e. scattering energies much lower than m^/M. At such 
energies the Compton wavelength of the nucleons is such that even the longest- 
ranged part of the NN force, that due to one-pion exchange, appears short-ranged, 
and so we can approximate the NN interaction by a potential of zero range. As we 
shall see this is an old story, dating back at least 50 years, and some of the recent 
work on EFT at these energies involved rediscovering the wisdom of the ancients. 
However, EFT provides new wrinkles to the old tales, and I will argue that these 
make the techniques developed by Hans Bethe and others all those years ago even 
more useful. 

2.1 Review of some basic NN facts 

To begin with, and to establish notation let us review some basic NN scattering 
facts. If the NN interaction V falls off sufficiently quickly at large r then the asymp- 
totic wave function has the form: 



C — Cd + Cd + I + Cd+2 + . . . + C-D+n- 



(44) 



2 The NN system in EFT(/r) 



^/>(r) = exp(ikiT) + /(A:,e) 



exp(ifcfr) 



(45) 



r 
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where kf = kf = k'^ = ME with E the energy in the centre-of-mass system, and 
cos^ = fci • fcf . The differential cross section for NN scattering is then given by 

^ = l/(fc,^)P. (46) 

Of course, if V is a central potential one may make a partial-wave expansion for /. 
In the absence of spin this has the form: 

f{k, 6) = ^ ^(2/ + 1) exp(i50 sinSi Pi{cos9) . (47) 

At this point the reader might reasonably complain that the NN problem involves 
spin and a non-central potential, but let me leave these complications aside for the 
moment, because for the ^Sq channel, which I will deal with first, neither actually 
matter. In this, or in any S-wave channel, the asymptotic radial wave function u 
has the form: 

uo{r) vo{r) ~ sin(fcr + 6o), (48) 

with 5q the S-wavc phase shift. The S-matrix itself is S = exp(2i5o) and this leads 
to a T-matrix (in the normalization used here): 

T = -i; , \ , . (49) 
M kcotS-ik ^ ^ 

This expression for the T-matrix incorporates the constraint of unitarity, i. e. 
probability conservation. An expansion for T, such as that obtained in EFTs like 
the one discussed above, only incorporates unitarity order-by-order in the EFT 
expansion. There is, however, a well-known expansion for the amplitude T which 
obeys (two-body) unitarity exactly. That is the effective range expansion (ERE), 
which is a power-series expansion for k cot 5. The ERE is a power series expansion 
for the real part of l/T. A simple derivation, valid for any energy-independent 
potential, shows that: 

k cot S{k) - lim k cot S{k) = ^k'^p{E,0), (50) 

where: 

1 f°° 

-p{Ei ^^2)= j [ve-, {r)vE2 (r) - ue-, {r)uE^ (r)] dr, (51) 
with UE{r) the S-wave radial wave function for the potential V at energy E, while 

sm()(fc) 

is a wave function that agrees with u in the asymptotic regime but differs from it 
inside the range of the potential, R. 
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Exercise 2 Derive Eq. ^5(\ ) by considering an energy-independent potential and 
analyzing the Schrddinger equation at two different energies Ei and E2. 

Thus, the region of support for the integration in Eq. (|5^) is restricted to r < R. 
This suggests that we may make a Taylor expansion for p in powers of E and that 
such an expansion should have a radius of convergence ^ 1 /R: 5 

p{E,0)^p{0,0)+E ^^f^^K .... (53) 
Combining Eqs. (|5^) and (^) yields Bethe's effective-range expansion ||ri|, [l^ : 

k cot 6{k) = -^ + ^rok^ - Pr^k'^ + ■■■■ (54) 

This can be converted to an expansion for 6 itself using the power series for the 
arccot, but the result only has a radius of convergence ~ 1/a. Thus, if i? <C a the 
effective range expansion is valid to much higher energies than a simple power-series 
expansion for 6. The effective range expansion was developed to treat precisely such 
a problem: NN scattering at low energies. 

Exercise 3 Show that for an energy-independent Hermitian potential of range R 
the effective-range is bounded from above by |Z^/; 

/ E>3 \ 

In the ""^So NN channel the central values of the effective-range parameters |^ 

are 

a (^So np) = -23.71 fm; 
ro (^So np) = 2.73 fm; 

P (^So np) = 0.06 fm, (56) 
while in the ■^Si channel we have jist : 

a f Si) = 5.42 fm; 

ro (3Si) = 1.75 fm. (57) 
The leading-order ERE amplitude is then: 



^ERE(o) _ 47ra 1 
Jnn - M 1 + iak' 



Meanwhile the second-order ERE amplitude is: 



ERE(2) _ _ 47r 1 



^) In fact the shape parameter P is not weU-determined by the experimental data, but it is 
definitely small. 
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Fig. 3. The effective-range expansion in the ^So channel (solid line) as compared to the 
Nijmegen phase-shift analysis (dots) [M. 



In the ^ So channel the amplitude ( |59| ) reproduces the experimental data up to mo- 
menta k ~ m^r (see Fig. even though strictly speaking the radius of convergence 
of the effective-range expansion is r7i^/2. 

The amplitude ( [s^ ) has a pole at momentum k = i/a. In the ^^Si system this pole 
corresponds to a low-energy bound state. This analytic structure of T(fc) must be 
reproduced if our effective theory is going to give a good account of NN scattering 
data in the low-energy regime. 

Exercise 4 Using the experimental numbers for the "^Si channel from Eq. ( |5^ 
compute the energy of the deuteron pole as given by the effective-range expansion. 
How does it compare with the experimental number: 

B = 2.224575(9) MeV? (60) 



2.2 Low-energy poles in EFT(yr) 

We now try to build an effective field theory to reproduce (^^ . Once this is achieved 
we will look at how to reproduce the NN amplitude if terms of higher-order in 
the effective-range expansion are incorporated, but we start with this very simple 
problem. The only degrees of freedom in our EFT will be the nucleons themselves, 
so we call this theory EFT(7f). The scale hierarchy in the theory is: 

k <C TO^, (61) 

where k is the nucleon momentum and l/m7r, which from now on will be denoted 
1/A, is the range of the force in the "underlying" theory. Note that here we have 
not initially specified whether the scales 1/a and l/rp are low-energy scales (like 
p) or high-energy scales (like A), although since 1/a « 8 MeV in the ^Sq channel, 
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it better ultimately be a low-energy scale! The symmetries of the theory will be 
translational invariance, rotational invariance (which may, however, be broken due 
to operators from one-pion exchange), and the Lorentz group for small Lorentz 
boosts. Only small boosts are allowed, since the theory is only valid for energies 
such that the nucleon momenta are much smaller than mjr, and so the nucleons are 
non-relativistic. Relativistic corrections can be incorporated systematically by de- 
manding that the theory is Lorentz covariant order- by-order in p/M. (We shall not 
discuss this further here, for details see Ref. (l|l.) With this in mind the Lagrangian 
for EFT(7f) is: 

C = N^ (^iD„ + ^'^N + Ci. (62) 

The interaction Lagrangian, Ci in principle contains all possible n-body interactions 
of nucleons, organized as a series in the dimensionality of the local operators that 
could contribute. For the two-nucleon system the lowest-dimension operator in Ci 
is (TV^iV)^, which has dimension six. Increasing numbers of derivatives and higher- 
body interactions can be included, but the corresponding operators are of higher 
dimension. Thus first we examine what we expect to be the "leading-order" (LO) 
theory: the NN scattering amplitude which arises from a Lagrange density with only 
a local momentum-independent four-nucleon operator. Therefore, the Lagrangian 
which we will use in our attempt to recover the amplitude ( p8| ) is: 

A = - [f^y^^'CoiN^PNy (N^PN), (63) 

where P is a spin and isospin projector that projects the NN pair onto the LSJT 
combination that is of interest to us — in this case, either the ^Sq or the "^Si, and 
n is the number of space-time dimensions. The factor (/i/2)^~" is included so that 
[Co] = —2 even when the theory is continued to fewer than four dimensions. 

Now, naive-dimensional analysis as formulated in the previous section implies 
that D^, as a dimension five operator, is suppressed relative to the dimension- four 
operator Dq. However, in a non-relativistic system the kinetic energy of our system 
is the same order as its total energy, and so D^/M is of the same order as Dq. 
Furthermore, in a non-relativistic bound (or quasi-bound) state such as we are 
considering here we also know that both of these quantities are of the same order 
as the potential energy. In this theory the NN potential is very simple. It is given 
by: 

(p'|y|p) = Co, (64) 

with p (p') the initial (final) relative momentum of the NN system. If we want 
V to be of the same order as the NN kinetic energy then the Co operator cannot 
be treated in the way suggested by naive dimensional analysis, since the naive 
dimensional analysis estimate is: 

Co - ^. (65) 
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x->o<->co<- ■ 

Fig. 4. Feynman diagrams contributing to the scattering amplitude from the Lagrange 
density in Eqs. (p2|) and (p3). Figure courtesy M. Savage. 



We would expect then, that Co, as a dimension-six operator, wiU be even more 
suppressed than . (The factor of AI appears here since this is a non-relativistic 
effective field theory 

This makes us wary of treating any contribution in the Lagrangian ( |63| ) in 
perturbation theory. But we will now see that the Lagrange density ( p3| ) is simple 
enough that the sum of all Feynman diagrams contributing to NN scattering, which 
is shown in Fig. ^, can be found exactly. 

Adding one additional loop to a Feynman graph adds one extra loop integration, 
one extra non-relativistic Green's function, and one extra Co vertex. Thus if k 
were the "typical" momentum scale of the problem, say the incoming NN relative 
momentum, then the higher-order Born series graph would indeed be higher order. 
It would contain an additional factor of 

n M 

fc^ X — X Co (66) 

as compared to its cousin which is one order lower in the Born scries for the potential 
(|6^. Putting in the NDA estimate ( p5| ) for Co indicates that adding a loop should 
reduce the size of the diagram by one factor of k/A. This is another way to see 
what we have guessed in the previous paragraph: if Co does scale according to 
Eq. (^5|) then there cannot be a low-energy NN bound state, since the potential 
( |6^ is "small" compared to the kinetic energy, and so the Born series can be used 
to calculate observables in perturbation theory. 

If, however, we blithely go ahead and sum the Born series depicted in Fig. ^we 
obtain a T-matrix 

Tnn — [ Co -I- CqIqCq + CqIoCqIqCq + . . . ] ~ —— — , (67) 

where Iq is the momentum-independent single-loop integral, which in the NN 
center-of-mass frame has the form: 

f d"(7 ( i 



2) J (27r)« \E + qo-(\^ /2Af + hi) \-qa-q^ /2M + v] 
/i\4-n r d("-i)q / 1 



2) J (27r)("-i) \E-(i^lM + iri 



(68) 



Here |k| — V ME the magnitude of the three-momentum of each nucleon in the 
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center-of-mass frame. Note that this integral Iq is linearly divergent, and so we have 
evaluated it in n space-time dimensions. 

In order to proceed further a regularization and renormalization scheme must 
be specified. It is clear from matching Eqs. ( |67| ) and ( |5^ ) to one another that the 
value of Co will depend on both the value of 1/a and the meaning assigned to the 
divergent part of Iq by our regularization and subtraction procedures. No matter 
how we define Iq, and hence the loop contributions, we can choose Cq so that ( |5^ ) 
is obeyed, but different regularization and renormalization procedures amount to 
a re-shuffling between contributions from the C'o vertices and the UV part of the 
loop integration (|6^). 

We begin by applying dimensional regularization with the MS subtraction scheme. 
Evaluating the integral in Eq. (|68| ) in the usual fashion of dimensionally-regulated 
loop integrals we see that 

lo = -M (-|kp - ie)("-3)/2 (^)'"" r (^^) (47r)(i-")/^ , (69) 

The integral has been regulated by the continuation to n dimensions. The linear 
divergence now manifests itself as a pole of the T function, which appears at n = 3. 
However, the vagaries of the F function mean that the n A limit of the expression 



in Eq. (69) is finite, so straightforward evaluation just gives: 



and thus the linear divergence in the integral never appears in the result for Iq. This 
is in contrast to the example in the previous section, where a logarithmic divergence 
in the integral in question revealed itself as a pole at the point n = 4 (e = 0) once 
the integral was continued to a different number of space-time dimensions. 

Equation (|7^) is the MS value for Iq. Re-examining Eq. (|67| ) we see that in this 
subtraction scheme each successive term in the bubble-sum involves the product 
CoM|k|. Thus, in order to produce a pole in the amplitude at momenta of order 
1/a it is necessary that the coefficient Co scale as Cq ~ a/M. In the "^Si channel 
this makes Cq significantly larger than its NDA value, and in the "'^Sq channel this 
makes Cq absolutely huge. This unnatural scaling in MS, and the problems it leads 
to, were first discussed by Kaplan, Savage, and Wise ||l9|. Cq is simply unnaturally 
large: in fact it is ill-defined for the case of a threshold bound state, since a ^ oo 
there. 

The real solution to this difficulty lies in the development of a new power count- 
ing |^0[ |2j, |lj, |2^, ^ . Subsequently it was shown how to implement this power 
counting on a diagram-by-diagram basis, using a general subtraction procedure 
called Power Divergence Subtraction |lj, ^ (PDS). In PDS both the n = 4 and 
n = 3 poles are subtracted from each bubble. The logarithmic divergence in n = 3 
corresponds to the power-law divergence in n = 4, and it then appears in the expres- 
sion for the bubble Iq. PDS thereby keeps track of linear divergences (see Ref. ||2^ 
for extensions of this idea) . This allows for "fine-tuning" between the coefficient Cq 
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and the linear divergence. In PDS, the integral Iq in Eq. ( |69| ) is defined to be 

a very similar expression to that obtained using a momentum cut-off, or with a 
momentum-space subtraction p9, GT\, ES 



If PDS is employed the bubble sum (p7[) yields 



Matching to the leading-order ERT amplitude (^8[): 

mCq a 

PDS is designed to isolate the linear divergence in the integral Iq. Ultimately 
this allows a portion of the attractive potential represented by Cq to cancel against 
the repulsive effect of the virtual kinetic energy of the nucleons in the potential 
well. The total effect of this kinetic energy is what is calculated in Iq. If the size of 
the well is R it scales as 1/{RM): (for the total number of states in the well) 
and 1/{MR?) (Heisenberg uncertainty principle) for the kinetic energy of particles 
in that well. By making Co, which really represents the depth of the well, also 
scale like R/M we can incorporate in our EFT the balance of potential and kinetic 
energy that we know is necessary to produce a shallow bound state. (See Ref. ||2^ 
for explicit implementations of this procedure using co-ordinate space regulation of 
the potential (p^).) PDS incorporates exactly this cancellation, but it also is more 
elegant than simply placing a co-ordinate space regulator at a distance i?, since the 
linear divergence is manifest within the framework of dimensional regularization. 

The explicit introduction of the subtraction scale into the expression for the 
amplitude results in Co having its scale set hy ^ 1/R) and not by 1 /a. Ultimately 
what this analysis shows is that in order for the four-nucleon interaction in Eq. (|63| ) 
to generate a system with a shallow bound state (i.e. 1/a ~ 0), Co must be near 
a non-trivial fixed point. The existence of this non-trivial fixed point 14 1 is in- 



dependent of the particular regularization and renormalization scheme chosen [p4[ , 
and its presence is what modifies the the scaling of Cq from the NDA ansatz (pq). 
Turning the argument of this section around we realize that i/we wish to construct 
a Lagrangian ( |63|) which reproduces the low-energy bound state that is present 
in the NN system then there must be a non-trivial fixed point in the flow of the 
coefficient in the LO interaction Lagrangian. 

Including higher-order terms in the EFT in the presence of such a fixed point is 
discussed in detail in Refs. |2^, ||. Here I will only say that to move beyond 

the leading-order ERE amplitude we design our EFT to reproduce the amplitude 
(E^ re-expanded as: 



T 



ERE(2) _ 



47r 1 



M-i + irofc2-ifc 
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Ana 
Ana 



1 + iak 
1 



M 1 + iak 



1 + iarofc^ 



1 



iak 
1 



1 + iak 



(74) 



In doing this we are assuming that ak ^ 1 but ar^k^ < |1 + iafc|. In fact this second 
condition is satisfied as long as k < l/rg. This technique works quite well in both 
the ^Si and ^Sq NN channels. Thus ultimately the scale hierarchy in the theory is: 



k ^ l/a = Q <t: A = l/ro 



(75) 



EFT(7f) expands low-energy NN observables in powers of Q/A, an expansion known 
as "Q-counting" . 

2.3 The deuteron in EFT(/) 

We will now see that a re-organization of EFT(7f) is useful if we wish to study 
deuteron properties. To do this we introduce an additional dibaryon field into the- 
ory, which we denote as t. The EFT(7f) Lagrangian (|6^) + (^3|) can then be rewritten 
as: 



ido 



2M 



N + it At - yit'iNN + tN^'N'']. 



(76) 



Note that the dimensions of the t field here are [t] — 3/2. This theory, or a close 
relative thereof, was proposed for NN scattering in Ref . |30j , but its antecedents are 
work due to Weinberg and the Lee model. The degrees of freedom are nucleons 
and the dibaryon field t. As in the Lee model there is are vertices for NN ^ t. 
Hence the dibaryon field must be dressed by NN loops (see Fig. ||). The Dyson 
equation for the dibaryon propagator Dt is then: 



1 

A 



-(-iSt)A 



with Et being able to be evaluated exactly in this theory: 



(77) 



(78) 



If MS is used to regulate Iq we discover that the t-field propagator for a dibaryon 
of energy E and momentum is just: 



DtiE;0) 



„,2 iMfc ' 

y 471- 



(79) 



with k^ = ME. Employing the relation: 

(-i)rNN = y^Dt (80) 

we find that the Tnn of Eq. ( ^ ) is recovered provided that the identification 



A 

y2 



n MS 



M 

Ana 



(81) 
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Fig. 5. The bare dibaryon propagator is dressed by nucleon bubbles to all orders, resulting 
in a dressed dibaryon propagator. Figure from Ref. S], courtesy S. Beane and M. Savage. 



is made. Note that this also shows that the theory ( |76| ) is equivalent to that (p^) of 
the previous subsection, provided the appropriate identification of A and y^vith 
Co is made. This can also be shown by constructing the path-integral for ( [76| ) and 
integrating over the field t to obtain the Lagrangian ( p2| ) + ( |63| ) p^ . 

(A technical note: it is, in fact, phenomenologically advantageous to adjust: 



My 



in the "^Si channel. Here 7 = y/MB is the binding momentum of deuterium. Nu- 
merically = 4.319 fm. By making this choice we ensure that the deuteron pole 
is at the correct energy.) 

A quick glance at the Lagrangian (|7^) reveals that the field t is not dynamic: 
it corresponds to a field of infinite mass, since it has no energy terms. If we extend 
the Lagrangian to include such a term: 



t-y[t^NN + tN^N^], (83) 



we observe that the field t actually has the wrong sign for its energy terms. This 
is a result of its auxiliary field nature. Repeating the derivation of Tnn from the 
previous paragraph we now obtain (using MS for the loop): 

rNN(i^;0)= ^_ . (84) 

My' Ait 

This gives us the opportunity to adjust the two independent parameters y and A 
to reproduce not just the first, but also the second, term in the effective range 
expansion for Tnn, i- e. to reproduce Eq. (|59|). 

Note that if we want to have both 1/a and ^rgk"^ in the denominator of our 
expression for Tnn then it must be that both are of the same order in our power 
counting. Otherwise we could re-expand the second-order effective-range amplitude 
as per Eq. ([z^). If we wish to have 1/a ~ ^rok^, and 1/a and k are both low-energy 
scales ~ Q then the logical conclusion is that 1/ro is also a low-energy scale ~ Q. 
In other words, both a and tq are unnaturally large and this should be taken into 
account in building the EFT ^ . 

Simple algebra then reveals that the choices: 

/ = A^s = (85) 
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will turn Eq. (p4) into Eq. (p9). Note that if 1/ro scales as Q, then y ^ Q 



1/2 



Exercise 5 Show that if PDS is used to regulate Iq the value of does not change, 
but: 



A 



PDS 



1 



- - 



(86) 



In retrospect it is better to write the Lagrangian for this theory in terms of a 
rescaled field t = yt. Then: 



ido + 



2M 



N - 



ido + 



4M 



- A 



[PNN + iN^N'']. 



(87) 



The dibaryon field i now scales with the small parameter Q as f ~ Q^/^Q^/^ = Q^. 

In the ^Si channel the dibaryon is a field for the deuteron, and the deuteron is 
"dressed" by NN loops, which simply means that unitarity is incorporated in the 
theory. In the case of the deuteron it makes more sense to rewrite Eq. (|59|) as 



T 



4^ 



-7^(fc), 



where, in the language of quantum field theory Zd is related to the wave-function 
renormalization of the deuteron field t. Meanwhile k ~ is the deuteron pole 
position in the complex fc-plane, and represents the "typical" momentum scale 
inside deuterium. 



(89) 



Exercise 6 Show that the second-order ERE amplitude ([5£) leads to 

-1/2 



1 



2ro 
a 



1.64. 



Hand-waving about power counting 

Effects due to higher-order terms in the ERE, mixing with the ^Di state, etc. 
can be systematically incorporated in the EFT (^). Here I focus on power counting 
for the coupling of the deuteron to photons. 

The single-nucleon coupling to photons is mandated by gauge symmetry and 
7N data. Meanwhile the unknown two-nucleon physics manifests itself in the theory 
as short-distance couplings of the photon to the dibaryon, for instance t^W'^Aot. 
Such effects are suppressed because "most" of the deuteron wave function "lives" at 
long-distances (see Fig. ^). Beginning with Eq. ( p8| ) and taking appropriate Fourier 
transforms one can derive the EFT(7f) deuteron radial S-wave wave function as: 

u{r) = y^2"fZd exp(-7r) . (90) 

This is, of course, exactly the same as the effective range theory (ERT) wave func- 
tion for this state. It stems from the idea that the long-range part of the deuteron 
wave function is prescribed by two numbers, the deuteron binding energy, which 
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Fig. 6. Radial S-wave and D-wave wave deuteron wave functions for EFT(7f) (or for ERT). 
Also shown is the wave function obtained from the Nijm 93 potential model [H. 



determines the exponential fall-off at long distances, and the coefficient of the expo- 
nential, which is known as in the literature and has an "experimental" value [ f^ 

As = 0.8845 ± 0.0008 fm~^/^ (91) 

This compares favourably with the ERT/EFT(7f) prediction from Eq. (^): As — 
0.871 fm"^/^ Effective range theory is an accurate way to extrapolate the NN 
amplitude from the physical region to the deuteron pole. 

In ERT, or in EFT(7f), this "asymptotic" wave function persists all the way in 
to r = 0, since in EFT(7f) the potentials involved are all zero range, and so nothing 
modifies the asymptotic form (|9^) until the origin is reached. We would expect that 
the asymptotic form would start to deviate from more detailed wave functions at 
r ^ 1/A — l/mTT, since that is the range of the underlying potential involved in 
NN physics. Thus the amount of wave function at short distances is parametrically 
small, being only 7/A of the total. However, in ERT there is no way to correct 
for this error, and so predictions of ERT should in principle be valid only up to 
corrections of relative size 7/A ~ 1/3. 

In the next subsection we will use the calculation of the charge form factor of 
deuterium to show that in EFT(7f) we can introduce operators which systematically 
correct for the (brutal) approximation that the form (|9^) persists all the way to 
r = 0. These operators are generically suppressed by (at least) one factor of the 
expansion parameter 

Q ^ (92) 

where q is the momentum of the external (electromagnetic or weak) probe of the 
deuteron. Q is thus the expansion parameter for calculations of processes on deu- 
terium in EFT(7f). For any given reaction the effects that can contribute are split 
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up into "long-distance" parts due to single nucleons interacting with an external 
probe at distances r 1/^,1/ q, and "short-distance" parts due to the probe in- 
teracting with the nucleons in the unknown region r ~ 1/A. The "long-distance" 
mechanisms can be calculated from knowledge of single-nucleon interactions with 
the external probes, and because of the "fluffy" nature of the deuteron wave func- 
tion (see Fig. |^) the short-distance operators are always suppressed relative to 
the long-distance ones. This gives EFT(7f) considerable predictive and explanatory 
power. 

2.4 Deuteron charge form factor 

The relevant Feynman rules for the calculation of the deuteron's charge form factor 
are as follows: 

Feynman rules and power counting: 

1. With each two-nucleon propagator we associate a factor: 



where E is the energy of the NN pair and p is the relative momentum. (This 
two-nucleon Green's function results from doing the relative-energy integra- 
tion for the product of two one-nucleon propagators.) Since E = —B ~ 
— both terms in the denominator scale with Q as /M , and so this 
propagator is assigned a factor Q^^. 

2. With each dibaryon-nucleon-nucleon vertex we associate a factor (— i). Clearly 
this carries no powers of the small momentum Q, and so scales as . 

3. Each loop corresponds to an integral over the loop three-momentum. (The 
zeroth-component integration has already been performed, see rule 1.) This 
results in a factor of . 

4. The interaction of an Aq photon with a single nucleon that results from 
the use of the minimal substitution Oq do — iQAq yields a jNN vertex 
corresponding to a factor — ie. We will count factors of e separately here, and 
this vertex then scales as e (surprise). 

5. The same minimal substitution in ( p3| ) yields a jii vertex ie/y^. Since y 
Q^^^ this coupling scales as e/Q and so is "super-leading", being enhanced 
over the coupling of a photon to a single nucleon. As will see this super-leading 
behaviour is a necessary consequence of electromagnetic current conservation. 

6. Other electromagnetic vertices not produced by the minimal substitution pro- 
cedure are also permitted in the theory. One such vertex involves a higher- 
derivative interaction of an Aq photon with a single-nucleon. This may be 
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represented as 

"finite nucleon size 



■^finite nucleon size = NW^AqN . (93) 



This generates a higher-dimension 'yNN vertex — ie(r^)q^/6. This scales as 
eQ^. (In fact, it is even more suppressed than this might indicate, since the 
isoscalar radius of the nucleon is not set by the pion mass, but instead by 
dynamics at the scale of chiral symmetry breaking A-^.) 

7. Analogously, there are vertices of the form e^i^V^Agi. Naive dimensional 
analysis, together with the knowledge that [i] = 2, suggests that [£] = —3. 
Hence we employ NDA to write ^ ~ 1/A^, and so have this vertex scaling as 
eQ^. The extra power of Q as compared to the scaling of the operator in (^) 
can be thought of as coming from the diffuse nature (technically known as 
"fluffiness" ) of the deuteron wave function. 

We are now in a position to calculate the charge form factor of the deuteron, Gc ■ 
In terms of standard quantum mechanics the charge form factor of the deuteron is 
defined by: 

Gc = f ((l|J°|l> + (Ok"|0) + (-lk"|-l)), (94) 



where we have labeled the (non-relativistic) deuteron states by the projection of the 
deuteron spin along the direction of the momentum transfer q and rj = |qp/ (4M^). 
Here it suffices to say that Gc(Q^) represents the amplitude for coupling of a 
"Coulomb" photon with four- momentum (0, q), = q^ to couple to the deuteron. 
Hence it is calculated by computing the amplitude for coupling of an Aq photon to 
the deuteron. 

The "leading-order" graphs for this process are depicted in Fig. ^ The graph 
on the right contains two two- nucleon propagators, one loop integral, and one 7NN 
vertex — ie. Hence it scales as: 

e 

e ^ —. 



Meanwhile the diagram on the left includes only the vertex ie/y'^ and so also scales 
as e/Q. Together they produce a deuteron form factor: 



sGc = y^7roZd 



(27r)3-B-p2/M' ' -B-{p + q/2y /M 

(95) 

where we have neglected deuteron recoil, and incorporated an overall factor of 
y'^jroZd in order to account for the wave-function renormalization that turns the 
amplitude for the AqU interaction into the amplitude for an Aq photon coupling 
to the deuteron. (For details on wave-function renormalization in this approach see 



Refs. |, 33|.) 
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Fig. 7. The leading-order diagrams that contribute to the deuteron charge form factor. 
The thick sohd hnes represent the deuteron field, while the thin lines with arrows denote 
the nucleon field. The crosses indicate the need for wave-function renormalization. Figure 
from Ref. [B3|, courtesy S. Beane and M. Savage. 



Exercise 7 Use the Fourier representation of 1/(MB + p^) in order to convert 
the integral in Eq. into an integral over r. You should find that 



167r2 



d?r 



exp(— 7r) 



exp(iq • r/2) 



exp(— 7r) 



- 1 



Note that the first term here can he interpreted as 

(V'|exp(iq- r/2)|V'), 



(96) 



(97) 



where li/") is the deuteron wave function corresponding to the radial wave function 
(^3). Evaluate the co-ordinate space integral to obtain: 



1 



"fro 



Aarctanf^l-l 

qra 



(98) 



The first correction to this result is due to the single-nucleon operator from 
Eq. (|9^). It gives an effect of 0{Q^) relative to those included in (p8|). This effect 
is easily included since the square of the nucleon's isoscalar charge radius, (r^), is 



known. The result 



is a prediction for Gc based only on low-energy NN data 



(ro and B) and current conservation. Current conservation imposes the constraint: 



Gc(g'-0)-l. 



(99) 



Note that without the leading-order effect of direct coupling of the Aq photon to 
the deuteron field this condition would have been violated, since the right-hand 
diagram of Fig. ^ yields for the deuteron charge, which is larger than one. This 
is just a consequence of the volume integral of the probability density corresponding 
to Eq. ( pO| ) being greater than one, but carefully maintaining current conservation 
in EFT(7f) ensures that Eq. ( p9| ) is satisfied. 

In EFT(7f) we obtain the 
(100) 



Meanwhile the slope of Gc at = gives ^^{r^) 



number: 



872 



(1.990 fm)2 



There are corrections to Gc, and hence to (r^), from short-distance operators in- 
volving Aq couphngs to tH, i. e. interactions of the last type discussed in the 
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Feynman rules listed above. These operators scale like g^, while leading-order pre- 
diction scales like 1/Q, thus we expect that the short-distance effects that are the 



first correction to (IOC) to be down by a relative factor of 7(7^ ^ Q'^. This means 



that the prediction (100) should be accurate to 3%. And, indeed this is the case 
with the difference between (IOC) and the experimental number (r^) = (1.971 fm)^. 
being on the order of 2% 



0. 



Although the result (100) can be derived in ERT too, systematic error estimates 
such as these are not possible in ERT. Also, it is not entirely clear how to include 
finite-nucleon-size effects in ERT, and the ERT wave function ( |90| ) will violate the 
current-conservation condition (p9[). As presented here, EFT(7f) is a descendant of 
Bethe's work, but it has evolved so that it has features that make it more suitable 
for the accurate computation of low-energy electromagnetic observables in the two- 
nucleon system. 

For what it is worth I have included a comparison of the EFT(7f) prediction 
discussed here, and computed in Refs. |3^, with the (sparse) low-g^ data on Gc 
(Fig. H). We see that the EFT(7f) calculation breaks down at momenta q ^ '^.m^^^ 
which is exactly what we would expect from the fact that at such momenta we 
begin to probe relative separations of the neutron and the proton where the theory 
completely fails to describe the dynamics In this regime we are no longer in the 
long- wavelength limit where the pion dynamics can be ignored. In Section ^ we will 
discuss how to improve the description of the deuteron state ji/)) in this important 
regime. 

While I do not have space here to do justice to the literature on EFT(7f) it 
is worth mentioning that a much more impressive and important application of 
EFT(7f') to low-energy electromagnetic processes involving the deuteron is the 1% 
calculation of the process np — > d7. EFT(7f) allows the derivation of analytic expres- 
sions for the El and Ml isovector strengths that dominate this process at low photon 
energies. Using the formalism discussed here, or the equivalent approach adopted 
by Rupak in his important calculation ||36|| we can derive results for 7d np that 
agree with the experimental data on low-energy deuteron photodisintegration at 
the expected level of accuracy. This is useful not only for what it teaches us about 
7d np, but also because it facilitates a computation of radiative neutron capture 
at the photon energies of relevance to big-bang nucleosynthesis by using a frame- 
work in which errors can be reliably estimated and analytic expressions developed. 
Similar work on vd interactions in EFT(7f') by Butler, Chen, and Kong, resulted in 
analytic expressions for the process v^d e~pp that were then used in computing 
the deuteron breakup rates critical to the assessment of the SNO data on neutrino 
oscillations |^ and also in calculating the solar hep process pp — > de^Vg |38j. Fi- 
nally, as a new generation of Compton scattering experiments on deuterium come 
on line EFT(7f) holds out the promise of empowering model-independent extractions 
of the neutron polarizability from the reaction 7d 7d p^. (See also below.) 



^) The factor of two in "g ~ 2m,r" arises because only half of the virtual photon's momentum is 
transferred to the deuteron's relative degree of freedom. The other half is deposited into motion 
of the deuteron's centre of mass. 
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Fig. 8. A comparison of Gc from EFT(7f) [||, H with a "fit" function (long-dashed hne) 
and experimental data points. The EFT(7f) calculation described above is represented by 
the dot-dashed line (i. e. it is NLO in the terminology of this figure), while the inclusion 
of finite-size effects for the nucleon then produces the, solid, "NNLO" curve, which is 
accurate up to corrections at relative order 0{eQ^) such as those due to short-distance 
effects in the deuteron wave function. Figure courtesy R. Gilman, from Ref. |35[ . 

To conclude this section, let me say that EFT(7f) as applied to processes on 
deuterium computes the "long-distance" portion of observables using wave func- 
tions that could have been written down using effective-range theory over 50 years 
ago. However, the errors in these wave functions can be systematically corrected 
for provided the probe is not of sufficient energy to truly "see" the "short-distance" 
part r < of the NN interaction. Thus, EFT(7f) is clearly a descendant of ERT, 

but it improves on it, since it allows for systematically- improvable calculations of 
(very) low-energy processes on the deuteron, and these resultant calculations have 
a well-justified theoretical error bar. This systematicity arises the existence of the 
small(ish) parameter 
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Such mingling of old and new to the benefit of nuclear physics is, I think, quite 
nice. To quote a very old source |4^, he (or she) who applies these ideas 

... is like the owner of a house who brings out of his storeroom new 
treasures as well as old. 

3 The three-body system with short-range interactions 

In this section we turn our attention to the three-body problem in EFT(7f). The 
problem to be solved is that of three particles interacting via forces that are short- 
ranged compared to the distance scales of interest. If these particles have spin-zero 
then one system where this physical condition is satisfied is the so-called helium 
trimer. Here the scattering length for two helium atoms is, 

02 ~ 125 1, (102) 

and these two helium atoms can come together to call a He- He bound state known as 
a dimer. Complicated Faddeev calculations suggest that the (three-body) scattering 
length for dimer-atom scattering is: 

03 ~ 195 A. (103) 

Both 02 and 03 are much larger than the typical range of the He-He interaction, 
which is about 10 A or so. 

As nuclear physicists we are also interested in the problem of neutron-deuteron 
scattering |^. Here the relevant scales are |l^, 42 : 



02 = 5.4194 ± 0.0020 fm; o^^^^^ = 6.45 ± 0.02 fm; o^^^^^ = 0.65 ± 0.04 fm, (104) 

where 02 is the "^Si scattering length. There are two three-body scattering lengths 
because the neutron is a spin-half particle, and the deuteron a spin-one particle, 
and so there are two possible channels for nd scattering: the S — 3/2 (quartet) and 
iS* = 1/2 (doublet), with S the total spin of the nd system. Since the range of the 
NN interaction is still R ^ I/ttItt ~ 1.5 fm, we hope and believe that EFT(7f) will 
be a viable way to treat 3N systems in low momentum reactions. 

As I will explain here, attacking the quartet channel in EFT(7f) turns out to yield 
remarkably accurate predictions for low-energy nd scattering without very much 
work. However, the doublet channel is much trickier, and significant conceptual 
work on the nature of renormalization in this problem was required before EFT(7f) 
could be made to yield answers. In fact, the conceptual problem in the 5 = 1/2 
channel also arises when discussing the scattering of three spinless bosons (such as 
helium atoms), and so, in order to avoid complications with spin, the helium-atom 
problem is the one we will discuss here. The good news is that this conceptual 

®) There is, of course, more proton-deuteron data, but in that case Coulombic effects somewhat 
mask the strong interaction effects we are trying to work on here. For a first effort at including 
Coulomb corrections in EFTs of the 3N system, see Ref. pc|]. 
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work has been (largely) done and there is now the promise of calculations akin 
to those described in the previous section, only for 3N systems. The discussion in 
this section mirrors closely that in the original papers on EFT(7f) applied to the 



three-body system. See especially Refs. 43, Bil Bsl M 



3.1 The equation 

We wish to describe the scattering of three spinless bosons via short-range interac- 
tions. Two of the bosons can form a bound state, which we include in the theory via 
an additional field t ( "dimer" ) . The Lagrangian for the theory is then exactly that 
of Eq. d?^ ) only with an additional term representing dimer-boson interactions: 

£ = iVt 1^0 + ^ N + t^At- y[t^NN + tN^N^]~ hi;^ijth. (105) 

(Note that due to a peculiar notational choice "iV" now represents a Helium atom. 
Note also that i has been discarded in favour of the original t.) We assume that 
all other higher-dimensional operators are indeed suppressed, and so work within 



the theory defined by Eq. (105). Furthermore, remembering that t is a dimension 
3/2 field, we realize that [h] = —2, so naive dimensional analysis implies that the 
three-body force term should be suppressed, and we thus expect to be able to set 
h — to get the leading-order result. We shall initially proceed in this way, but it 
is worth warning that we have already seen that in the presence of "unnaturally" 
long distance scales, such as the scattering length ( |102D , NDA is too naive. 

Now consider boson-dimer scattering within this theory. The incoming dimer 
is the dressed field, with the propagator given by Eq. (^9|). All of the graphs that 
contribute to boson-dimer scattering in the theory with h = are shown on the top 
line of Fig. ^ These graphs just form the multiple-scattering series for the three- 
body problem with short-range interactions, and can be re-summed by writing the 
boson-dimer scattering amplitude as the solution of a Faddeev equation. Indeed, 
in this case, since the potential is separable the Faddeev equation reduces to a 
Lippmann-Schwinger equation with one-nucleon exchange as the potential (bottom 
line of Fig. I). 

This can be understood from an EFT point of view by simple power counting. 
Each additional loop in one of the graphs on the top line of Fig. ^ adds two tNN 
vertices (each ~ y), one extra non-relativistic three-nucleon propagator (^ M/Q'^), 
and one extra tN propagator (~ l/{My'^Q)). Including the three-dimensional loop 
integration we see that the overall effect is that a loop is down by: 

i. e. it is not suppressed at all. Thus it is necessary to sum this entire class of loop 
graphs just to get the LO answer for the boson-dimer scattering amplitude. 

To derive this equation we choose the kinematics shown in Fig. ^j. The outgoing 
(incoming) dimer field carries momentum k (p). Meanwhile, the outgoing (incom- 
ing) boson carries momentum — k (— p), The outgoing deuteron and nucleon are on 
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Fig. 9. The first line shows the infinite sum of graphs giving the leading order neutron- 
deuteron scattering in the spin 3/2 channels. The second line is the graphical representa- 
tion of the Faddeev equation that sums them up. Figure from Ref. 



shell, and e denotes how far the incoming particles arc off shell. This is the "half- 
off-shell" amplitude for tN scattering. We now let i(p, 0) be the amplitude for the 
sum of all Feynman graphs in the theory (|7^) which contribute to tN scattering in 
this kinematics (with h = here): 

The integral equation for t can then be read off from Fig. ^. At k = it is: 

'(P.O) = ^.«^a/^ . . (107, 

where the artificial parameter A has been introduced, with A = f for the case of 
boson-dimer scattering. 

Exercise 8 Check that applying the Feynman rules for the field theory nm) leads 



p2 + 72 J (2^)3 72 + p2 + q2 + p . q 



to Eq. I\101^ . 

At zero relative momentum only S-waves contribute and so we postulate: an 
amplitude 

/ „x 3Af 1 , , M , , , , 

t(p,0) = -- / a{p) = ——a(p). 108 



Working with a{p) rather than 0) makes our life easier when we come to compare 
with experiment, since a(0) = —03, the negative of the three-body scattering length. 
Equation (107) now boils down to a one-dimensional integral equation for a{p): 

■ a[p) = ^- — J + — log / dq — log „ „ — 5 a{q) . (109) 



F(j>) p^ + TT J 2pq \p'^ + q"^ + — pq ^ 
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This is the Faddeev equation for S-wave scattering for the case of contact forces 
derived firstly in reference |M by different methods. (Note: in fact, it was actually 



derived before Faddeev derived his equations!) Eq. (109) can now easily be solved 
numerically. 



Exercise 9 Derive Eq. (UjMI from Eq. (101). 



3.2 The problem 

Our next goal is to derive the ultraviolet (i.e. high momentum) behaviour of the 
half- off-s hell amplitude a(p). For p S> 7 the second term on the right-hand side of 



Eq. ( 109 ) dominates, and the integral equation becomes: 



a{p) - -j=- / — log — — ^ . 110 

A few things about this equation are worth noting: 

1. A cutoff has been applied, in order to regulate otherwise troublesome be- 
haviour from the high-momentum part of the integral. 

2. This is a homogeneous integral equation. No information from the driving 
(inhomogeneous) term survives in the UV limit. Some authors have argued 



that this makes the original problem defined by the equation (109) ill-posed, 
since the kernel is non-compact [^] , but in fact the kernel is compact provided 
that any finite cutoff is imposed. 



3. If a{p) is a solution of Eq. (110) then so is its complex conjugate a*{p). 



4. If we were to take A — > 00 then the integral on the right-hand side is scale 
invariant. 

5. Similarly, if A — > cxd and a(p) is a solution then a(l/p) is also a solution. 

This last point suggests that we try power-law solutions of the form a{p) — . 
Doing this, and performing the resulting McUin transform on the right-hand side, 
we reproduce an old result of Danilov |^ : The homogeneous equation ( |llO| ) will 
always have such a solution provided that s obeys the transcendental equation: 

^_ 8^sin(W6) 0^ (111) 
y/3s cos(s7r/2) 



Now, the solution a{p) has very different properties depending on whether Eq. (Ill) 
admits real or complex solutions for s. Thus, we divide our discussion into two cases: 
A < Ac = in which case the solution of (111) is real, and A > Ac, in which case 



it is complex. Note that the case we were originally considering was A = 1 > Ac 
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3.2.1 The simple case: A < Ac 

However, t he ca se A < Ac is also interesting, not least because it transpires that 
the equation ( |109| ) with A = —1/2 describes neutron-deuteron scattering in the 
quartet (spin-3/2) channel. In this case solutions to Eq. (] 



111 



are 



s = ±2,±2.17, .... (112) 

Here we discard the solutions that grow with p as we go far off-shell, and so obtain: 

a{p) = Cp-\ (113) 

It follows that the integral equation is well-behaved and the results are completely 
insensitive to the size of the cutoff A, and to the form of the cutoff function em- 
ployed. (As long, of course, as A » 7.) We could express this schematically by 
saying that: 

da(p) 



A 



dA 



«0. (114) 
low p 



This allows us to make a confident prediction for a{p), and hence for 03 = — a(0). 
That prediction is that the quartet neutron-deuteron scattering length is [pj[: 



3/2 _ 



5.09 fm. (115) 



However, we should remember that the deuteron propagator employed to obtain 
this number reproduced only the deuteron binding energy. Putting in a deuteron 
kinetic energy, allows us to also reproduce As, or, equivalently, tq. If this is done, 
and the field theory (p3) used to perform the above analysis we instead find [p3[: 



3/2 



6.33 ± 0.1 fm; (116) 



(See Ref. |50) for an earlier, similar approach to this problem.) Note that the cor 
rection due to these "range effects" is about 20% of the leading-order 03^^ ( 115 ) 



a number that might have been anticipated given the fact that "fr^ ~ 0.4. The 
error ±0.1 fm comes from using analogous scale arguments to estimate the effect 
of including more information about NN scattering in the deuteron propagator, as 
well as what would happen if NNN forces were to be included in this channel. 



The number (HE) compares well with the experimental value and with sophis- 



ticated Faddeev calculations using modern, accurate, NN potentials M. In fact. 



the calculation (116) is model independent, and involves only the input parameters 



7 and ro, thus one would expect any potential model which reproduces these num- 



bers to agree with the result (116), within the estimated error. It is also important 
to note that model-independent statements of similar accuracy can be made for 
the energy-dependence of the quartet nd S-wave phase shift |Q and for higher 
partial- waves (in both the quartet and doublet channels actually) [ p2[ . 
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Fig. 10. The amplitude a{p) (denoted here K(0,p)), in the bosonic case for different 
values of the momentum cutoff A in units of the two-body scattering length 02. Figure 
from Ref. K^, courtesy U. van Kolck. 



3.2.2 The tricky case: A > Ac 

The "problem" alluded to in the title of this subsection does not arise until we 



consider couplings A sufficiently large to generate imaginary solutions to Eq. ( |lllD . 
In particular, for the case originally of interest here, A = 1 the solution is: 

s = ±iso, (117) 

with So — 1.0064. Then, in the intermediate domain between the low- momentum 
scale 7 and the cutoff A we anticipate that the solution a{p) will behave as: 



a(p) = Ccos[solog(^) +5] . (118) 



Since there are now two admissible solutions to the integral equation (|110|) two 
parameters which are not determined by the integral equation make an appearance 
in a{p): C and 8. The value of these parameters is ultimately fixed by dynamics 
near the cutoff A. And this means that the phase 5 is fixed by ultraviolet, or short- 
distance, physics. The result is that in this case the low-momentum behaviour of 



a{p) varies wildly as the cutoff is changed. (See Fig. 10.). This is unacceptable from 
a theoretical point of view: the cutoff is a theoretical artifice, and dynamics near 
the cutoff should not play a key role in determining the low-energy scattering in 
general and the scattering length in particular. 

One way to understand this difficulty is to realize that the nucleon-exchange 
potential here can be thought of as a 1/r^ force law if we work in hyperspherical 
co-ordinates and consider 1/A <C r <C I/7. This is a singular potential, and the 
corresponding Hamiltonian is unbounded from below. The difficulties in defining 
quantum-mechanical problems based on such interactions are well-known |5^, p4[ , 
since systems tend to "fall to the centre" as they can always lower their energy by 
cramming more of the wave function into the short-distance region. Of course, this 
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Fig. 11. Equation for the three-body amplitude including the three-body force. Figure 

from Ref. [Bl. 



is not advantageous if we are dealing with a repulsive l/r"^ potential, and so that 
problem is well-defined. Now we gain some insight into why the nd quartet resulted 
in a well-defined quantum mechanics problem. In that case the Pauli principle 
precludes all three nucleons from getting close together, and the result is a repulsive 
nucleon-exchange potential between the deuteron and the third nuclcon. In contrast, 
in the triton, or in the three-boson system, the Pauli principle does not forbid 
"collapse to the centre" , and both the triton and the trimer cannot be described by 



Eq. ( 109 ), since that equation fails to make a concrete prediction for 03, or indeed, 



for any low-energy scattering observable. 



3.3 The solution 



The presence of such sensitivity to short-distance physics suggests that the problem 
requires renormalization. For those used to solving EFTs in perturbation theory 
this is surprising, since the difficulty the Eq. ( |109| ) does not arise at any finite or- 
der in perturbation theory, it only occurs when we consider the non-perturbative 
solution to the problem of neutron-deuteron scattering Efl, Esl. However, in con- 



sidering the non-perturbative solution to Eq. (109) it is perhaps not surprising, 
given non-compact nature of the integral equation's kernel in the limit A — s- cx), 
or, equivalently, the singular behaviour of the one- nucleon-exchange "potential". 
Regardless of whether it is surprising or not, the presence of the unphysical be- 
haviour with variation in A depicted in Fig. |l^ suggests that we add to the kernel a 
counterterm. A natural candidate is the three-body interaction we already included 
in the Lagrangian ( |105| ). In retrospect, dropping the three-body force term from 
the integral equation is not justified in the case of bosons (although it is justified 
in S-wave quartet nd scattering). Graphs containing the three body force are now 
included as shown in Fig. (pi]). 

The new graphs depicted here modify the kernel of the integral equation for 
threshold scattering to: 

K{q,p) ^ — \og\— , 2 I 2 + T2' 



where H is related to the coefficient h appearing in the Lagrangian ( 105 ) by: 
so that the parameter H is dimensionless. 



(120) 



36 
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Fig. 12. The running of t he " three-body force" H{A) as a function of A. The sohd hne 
is the prediction of Eq. (123) and the points are values obtained for H by adjusting 
it numericaUy so as to reproduce a given three-body scattering length as. Figure from 

Ref. Hi courtesy of U. van Kolck. 



For H{A) ^ 1 the three-body force has an effect in Eq. (11£) only for p in 
the vicinity of the cutoff. Therefore the analysis of the behaviour of the amplitude 
given in the previous subsection is still valid in the intermediate region 7 ^ p ^ A. 
It is now relatively straightforward to choose H = H{A) in such a way that the 
low-momentum part of a{p) is (essentially) independent of cutoff. In other words 
we want to try and effect a choice of H such that: 



a{p) = C{A, H{A)) cos [so log (^) + S{A, H{Aj)] 



(121) 



contains a phase which is independent of the cutoff. To this end we define an 
equation for the evolution of H{A) with the cutoff A (a "renormalization group" 
equation) : 

- so log A + 5{A, H{A)) = -so log A, , (122) 

where A^ is some fixed physical scale. 

A fairly simple calculation (see Ref. |4^) then reveals that an (approximate) 
solution to this is provided by: 



HiA) = 



sin (so log [A/ A*] — arctan [1/so]) 
sin (so log [A/A,] + arctan [1/so]) 



(123) 



Choosing H as per this formula, then inserting it into the kernel (119) and solving 
the integral equation for a{p) should ensure that the low- momentum behaviour of 
a is relatively insensitive to the cutoff A. This procedure can also be implemented 
numerically, and the results agree startlingly well. (See Figure p^ ) 

The key point here is the appearance of a new physical parameter (e. g. A, 
here) that is not determined by two-body physics. This means that in systems such 
as the Helium trimer, or the nd doublet channel, on-shell two-body scattering data 
alone is not enough to fix the three-body observables. One additional parameter 
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Fig. 13. Predictions for k cot S as a function of the relative momentum of the dimer and 
heUum atom k, in the case that as = 1.52a2. The different curves are for four different 
cutoffs (Aa2 = 42.6,100.0,230.0,959.0), with H adjusted in each case to ensure that 
as = 1.52a2. Figure from Ref. courtesy of U. van Kolck. 



is required, and this can be thought of as A* itself or the value of the three-body 
force iJ(A*) at that scale. 

Here I have not discussed the integral equation for the half-ofF-shell scattering 
amplitude at finite energies. Sufficed to say that this can be derived in exactly the 
same fashion as was done for the zero-energy equation here. And the issue of the 
low-momentum amplitude's sensitivity to the cutoff also appears at any finite en- 
ergy. However, in the region /c ^ p <C A (fc the on-shell momentum of the dimer 
and boson) the above analysis applied to a{p) goes through, since all reference to 
the scale k drops out. Consequently, the renormalization of the problem at k — 
through the addition of an energy-independent H(A) is sufficient to renormalize 
the problem at any finite k too. (For a nice way of viewing this see Ref. [^ .) Thus 
the energy-dependence of the phase shifts is predictive. In boson-dime r sca tter- 
i ng r cnormalizing in order to reproduce the ratio ajijai given by Eqs. (|l03| ) and 



(102), yields the curve shown in Fig. |l^ for A; cot (5 in S-waves. The curves shown 
are for different choices of the cutoff, with i/(A) adjusted in each case to give the 
same A;cot(5|j,^g. Once this has been done remaining effects of the regulator are 
suppressed by k^ j ts? . These results are actually in accord with the (scarce) experi- 
mental data on the trimer system, and also with sophisticated model calculations of 
Helium atom-dimer interactions. The beauty of the EFT approach described here 
lies in its simplicity, and its consequent ability to generate universal predictions for 
energy-dependence of kcoXb once the ratio a^ja-i is specified. 

Such a universal prediction can be made for the relationship between any two 
low-energy observables. For instance, once 03 is fixed, the three-body bound-state 
energy B3 is a prediction of the theory. For a given choice (A, 77(A)) we have a 
definite prediction for both 03 and ^3 and so we can generate a parametric plot of 
the function ^3(03). This is what has been done in Fig. |l^. In this case the plot is for 
the case of the triton, and I will suppress the additional spin-isospin complications 
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Fig. 14. Results for the doublet S-wave neutron-deuteron scattering length for different 
models having very similar two-body scattering lengths and effective ranges, denoted by 
the dots. The solid line is the EFT prediction. The cross indicates the experimental value. 
Figure courtesy U. van Kolck, from Ref. [p6|. 



that arise in this channel. (For details see Ref. |Q .) I will just point out that exactly 
the physics I have been discussing here results in the one-parameter relationship 
between and 03. The existence of such a relationship is an old story in the nd 
system Q Q]. Plotting the predictions of a number of different NN models which 
all give the same low-energy NN data, we see that they give different predictions 
for both i?3 and 03. But these different predictions are absolutely correlated: if we 
put all of them on a graph of i33(a3) then they fall on a straight line, known as the 
Phillips line (not me, another Phillips). 



The plot Fig. 14 shows that the existence of this Phillips line can be understood 



as being due to the dual facts: 

— Two-body NN data is not sufficient to determine the low-energy scattering in 
the doublet nd channel. 



— One, and only one, additional parameter, which can be adjusted to one piece 
of NNN experimental data, then allows concrete predictions to be made for 
low-energy NNN observables. 



^) The Phillips line has traditionally been understood as arising from different nff-shell be- 
haviours in the two-body NN amplitudes represented by the different NN models pq . In fact, 
this explanation is not in conflict with the one given here, in which the Phillips line is under- 
stood through the necessity for a three-body force, since field redefinitions change both off-shell 
behaviour of amplitudes and strengths of three-body forces |59[ . 
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4 Adding pions: a brief discussion of using chiral perturbation theory 

for NN-system processes 

We now wish to discuss energies larger than the very low energies at which 
the effective field theory without explicit pion degrees of freedom is valid. We will 
consider interactions of nucleons and pions in the kinematic regime where 

p q ^ ^ nip, M (124) 

with p the (three)-momentum of the nucleons, q the (four)-momcntum of any ex- 
ternal probes, nip the rho-mcson mass and AI the nuclcon mass. 

In this kinematic domain an important constraint on pions' interactions with 
each other and with nucleons is provided by the spontaneously-broken approximate 
chiral symmetry of QCD. I will not provide a review of the physics of chiral sym- 
metry here. Sufficed to say that extensions of the techniques described in Section 
1 can be used to derive an effective field theory in which the low-energy degrees of 
freedom are nucleons and pions. The symmetry which constrains their interactions 
is the approximate S'[/(2)l x SU{2)y{ symmetry of QCD. But since this symmetry is 
spontaneously broken and approximate the development of this effective field theory 
requires that we also understand how to introduce into the EFT operators which 
break 5f7(2)L x SU{2)jx in the same fashion in which it is broken in QCD. The re- 
sultant Lagrangians, the development of which is reviewed in Refs. and 
in the lectures of Prof. Leutwyler in this volume, represent an effective quantum 



field theory valid in the kinematic domain (124). The resulting theory is known as 



chiral perturbation theory (xPT). The small parameter in which we expand is now 

^-A^'A^'f^' 

A^SB A^SB JVjj-sb 

where A^sb is the scale of chiral symmetry breaking: 

A^SB^mp,M,47TU. (126) 

To summarize years of work on xPT in a (very) hand-waving fashion, the two 
key ideas that make the xPT expansion tenable are: 

— That the short-distance effects of QCD which are not strongly constrained 
by chiral symmetry can be expanded in powers of the small momenta in 
the problem p and q, with the resulting operators having coefficients whose 
sizes are determined by naive dimensional analysis, with the high-energy scale 
being A^sb- 

— That the "long-distance" loop effects due to low-energy pions and nucleons 
interacting involve only weak interactions, since chiral symmetry implies that 
the chiral Lagrangian can be organized in such a way that the pions only 
ever couple derivatively to the nucleons and to each other. Thus every pion 
interaction involves a factor of 

^, (127) 
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and so pionic loops are ultimately suppressed by factors of 

{k,m^) 



(128) 



where k is the magnitude of the pion's three momentum. 

We have already seen that the "typical" momentum scale relevant for binding 
in deuterium is \/ MB which is small compared to A^ge- Although the typical 
momentum in the trinucleons is higher, it is still well below A^sb- Thus, we expect 
to be able to calculate the response of such nuclei to low-energy probes using this 
effective field theory. The result is a systematically-improvable, model-independent 
description of electroweak processes on these nuclei. Furthermore, pionic processes 
now fall within the purview of the EFT, and although I will not discuss them here, 
calculations of reactions such as 7rd nd and 7d n'^d in xPT have been quite 
successful ^ ^ . 

In this section I will discuss two recent calculations of electromagnetic reac- 
tions on deuterium in chiral perturbation theory that I have been involved in. The 
reactions are elastic electron-deuteron scattering and Compton scattering on deu- 
terium. 

4.1 Power counting 

We begin by establishing the power counting for processes in this approach. While 
this power counting can be expressed diagrammatically, a la that derived above for 
EFT(7f), here I will lay out the power counting in the fashion of a simple quantum 
mechanics problem. If you want to derive the diagrammatic rules from the ones I 
give here I leave that to you as an exercise! 

Consider an elastic scattering process on the deuteron whose amplitude we wish 
to compute. If O is the transition operator for this process then the amplitude in 
question is simply (7/;|O|'0), with \^p) the deuteron wave function. Here I follow 
Weinberg j6^, ^ , and divide the formulation of a systematic expansion for this 
amplitude into two parts: the expansion for O, and the construction of lip). 

Chiral perturbation theory gives a systematic expansion for O of the form: 



6 = ^0("), (129) 



n=0 

where we have labeled the contributions to O by their order n in the small parameter 



P of Eq. (125). Eq. ( |129| ) is an operator statement, and the nucleon momentum 
operator p appears on the right-hand side. However, the only quantities which 
ultimately affect observables are expectation values such as {tp\p\ip). For light nuclei 
this number is generically small compared to A-^sb. 

To construct O^"-* one first writes down the vertices appearing in the chiral La- 
grangian up to order n. One then draws all of the two-body, two-nucleon-irreducible, 
Feynman graphs for the process of interest which are of chiral order P". The rules 
for calculating the chiral order of a particular graph are: 
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— Each single- nucleon propagator scales like 1 /P, since in (heavy-baryon) xPT 
the nucleon is a static object (at leading order) and its propagator is just 
1/po, with pq the amount by which the nucleon's four momentum is taken 
away from the mass-shell by the process in question; 

— Each loop contributes P**; 

— Graphs in which both particles participate in the reaction acquire a factor of 

— Each pion propagator scales like 

— Each vertex from the nth-order piece of the chiral Lagrangian contributes P". 

These rules are consequences of naive dimensional analysis, pure and simple. The 
argument is simply that the powers of P derived from considering the construction 
of 6(") in this way must be carried by the low-energy scales p, q, and m^. Meanwhile 
as the dimension of O*-") increases more and more powers of the high-energy scale 



Ajj-SB must appear in the denominator in order for Eq. (12£) to be dimensionally 
correct. 

While this is a simple argument it has the significant consequence that more 
complicated graphs contributing to O, i. e. ones that involve two-body mechanisms, 
and/or higher-order vertices, and/or more loops, are suppressed by powers of P. 

4.2 Deuteron wave functions 

There remains the problem of constructing a deuteron wave function which is con- 
sistent with the operator O. Weinberg's proposal was to construct a xPT expansion 



in Eq. (129) for the NN potential V, and then solve the Schrodinger equation to 
find the deuteron (or other nuclear) wave function ^ Recent calculations 
have shown that nuclcon-nuclcon scattering data can be understood, and deuteron 
bound-state static properties reliably computed, with wave functions derived from 



xPT in this way |69|, |70| [n] [72[ [73|. 

Now for xPT in the Goldstone-boson and single-nucleon sector loop effects 
are generically suppressed by powers of the small parameter P. This is just a 
consequence of the weakness of pionic interactions in the chiral limit. Thus, in 
these reactions the power counting in P applies to the amplitude, and not to the 
two-particle potential. However, the existence of the deuteron tells us immediately 
that a power counting in which loop effects are suppressed cannot be correct for 
the two-nucleon case, since if it were we could use successive Born approximations 
to compute the NN scattering amplitude, and there would be no NN bound state. 
Weinberg's proposal to instead power-count the potential is one response to this 
dilemma. However, its consistency has been vigorously debated in the literature 
(see [0, H for reviews). Recently Beane et al. Q have brought some clarity to 
this discussion, by showing that Weinberg's proposal is consistent in the ^Si — ^Di 
channel. In fact, pionic interactions are not "weak" in this channel since even in 
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Fig. 15. Radial S-wave and D-wave wave deuteron wave functions for several different 
NN potentials. The dotted line is the ERT/EFT(7f) wave function, while the dot-dashed 
line is that from a simple model with one-pion exchange at long distances and a square 
well at short distances Q. The long-dashed line is the NLO deuteron wave function of 
Ref. and the solid line is the Nijm93 wave function [M. 



the chiral limit a (tensor) potential 

V^tW ~ ^, (130) 

exists between the two nucleons. This is an (infinitely!) strong interaction. Beane 
et al. propose that xPT be reformulated in NN systems as an expansion about the 
chiral limit. They have shown that this is a consistent proposal, and that there are 
some signs that such an expansion will converge, albeit not particularly rapidly. 

Here I adopt a more "lowbrow" approach. One way to understand xPT power- 
counting for deuteron wave functions is to examine the deuteron wave function in 
three different regions. 

1. In the region R 3> I/tItt the deuteron wave function is described solely by 
the asymptotic normalizations Ag, A]j, and the binding energy B. These 
quantities are observables, in the sense that they can be extracted from phase 
shifts by an analytic continuation to the deuteron pole. This region is the one 
in which EFT(7f) gives the "correct" deuteron wave function. 
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2. The second region corresponds to R ^ l/rriTr. Here pion exchanges play a 
key role in determining the NN potential V , and, associatedly, the deuteron 
wave functions u and w. The leading effect comes from iterated one-pion 
exchange — as has been known for at least fifty years. Calculations with onc- 
pion exchange (OPE) defining the potential in this regime will be referred 
to below as "leading-order" (LO) calculations for the deuteron wave func- 
tion. Corrections at these distances come from two-pion exchange, and these 
corrections can be consistently calculated in xPT. Two-pion exchange effects 
are suppressed by powers of the small parameter P, with the "leading" two- 
pion exchange suppressed by relative to OPE. This two-pion exchange 
can be calculated from vertices in C^^^ and its inclusion in the NN potential 
results in the so-called "NLO" calculation described in detail in Ref. [ [7l| . 
Corrections to this two-pion-exchange result from replacing one of the NLO 
two-pion-exchange vertices by a vertex from C)^-^- This results in an additional 
suppression factor of P, or an overall suppression of P^ relative to OPE, and 
an "NNLO" chiral potential g |7^, |7l|, |7|. 

3. Finally, at short distances, R <C '^/m-^ we cutoff the chiral one and two-pion- 
exchange potentials and put in some short-distance potential whose param- 
eters are arranged so as to give the correct deuteron asymptotic properties. 
This cutoff will generally be kept finite here, ft remains an open question 
whether it is useful to attempt to take i? — > 0, in the manner that we essen- 
tially did for EFT(7f). 

A plot of different deuteron wave functions is displayed in Fig. 



4.3 Elastic electron scattering on deuterium 

Both the deuteron's charge form factor, Gc (defined in Eq. (^)), and its quadrupole 
form factor, Gq, involve matrix elements of the zeroth-component of the deuteron 
four-current, J*^. Here we split into two pieces: a one-body part, and a two-body 
part. The one-body part of begins at order e, with the impulse approxima- 
tion diagram calculated with the non-relativistic single-nucleon charge operator 
for strutcutreless nucleons. Corrections to the single-nucleon charge operator from 
relativistic effects and nucleon structure are suppressed by two powers of P, and 
thus arise at O(eP^), which is the next-to-leading order (NLO) for Gc and Gq. 
At this order one might also expect meson-exchange current (MEC) contributions, 



such as those shown in Fig. 16. However, all MECs constructed with vertices from 
'^ttJV isovector, and so the first effect does not occur until N^LO, or O(eP^), 

(2) 

where an NNttj vertex from C^^p^ replaces the upper vertex in the middle graph 
of Fig. and produces an isoscalar contribution to the deuteron charge operator. 
(This exchange-charge contribution was first derived by Riska [75[|.) 

The most important correction that arises at NLO [0(eP^)] is the inclusion of 
nucleon structure in xPT. At this order isoscalar form factors are dominated by 
short-distance physics, and so the only correction to the point-like leading-order 
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Fig. 16. The impulse-approximation contribution to Gc and Gq is shown on the left, 
while two meson-exchange current mechanisms which would contribute were the deuteron 
not an isoscalar target are depicted in the middle and on the right. 

result comes from the inclusion of the nucleon's electric radius, i.e. 

4^^XPTNL0 = 1-^(4^^V. (131) 

This description of nucleon structure breaks down at momentum transfers q of order 
300 MeV. There is a concomitant failure in the description of ed scattering data (see 
Figure |l^). (This failure was already seen in the work of Meifiner and Walzl |^.). 
In order to focus on deuteron structure, in the results presented below I have chosen 
to circumvent this issue by using a "factorized" inclusion of nucleon structure Jt?] . 
This facilitates the inclusion of experimentally-measured single-nucleon form factors 
in the calculation, thereby allowing us to test how far the theory is able to describe 
the two-body dynamics that takes place in ed scattering. 

We now compute results for Gc and Gq at LO (structureless nucleons), NLO 
(adding nucleon structure and the first relativistic corrections), and N^LO (includ- 
ing the meson-exchange contribution discussed above in the charge operator Jq). 
This yields the parameter-free predictions shown in Fig. |l8|. The figure demon- 
strates that convergence is quite good below q ~ 600 MeV — especially for Gc ■ The 
results shown are for the NLO chiral wave function, but the use of the NNLO chiral 
wave function, or indeed of simple wave functions which include only one-pion ex- 
change, do not modify the picture greatly below q — 700 MeV [Q. It is also clear 
that — provided the experimental data from electron-nucleon scattering is taken into 
account — xPT is capable of describing the charge and quadrupole form factors of 
deuterium at least as far as the minimum in Gc ■ This result is encouraging for the 
application of EFT to light nuclei. 



4.4 Compton scattering on deuterium 

The scattering of low-energy photons from the nucleon at low energies probes the 
long-distance structure of the neutron and the proton. In the case of the neutron 
the low-energy Hamiltonian for the neutron's interactions with photons is exactly 
(|l|). The coefficients of the and tcr ms are related to the electric (an) and 
magnetic (/3n) polarizability of the neutron via Q|: 

ai — — 27ran; 02 = — 27r/3n. (132) 
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q (MeV) 

Fig. 17. The deuteron charge form factor at orders e, eP^ and eP'^ in a "strict" chiral 
perturbation theory calculation. The experimental data is taken from the compilation of 

Ref. 



Photon- nucleon scattering has been studied in xPT in Ref. ||79|, |80[] , where the 
following results for these neutron polarizabilities were obtained at leading non- 
vanishing order (O(e^P)): 



- 12.2 X 10-4 fm^; /3„ = ^Jff; = 1.2 x 10"^ fm^ . (133) 



Chiral perturbation theory also predicts an = Op, /?n = /3p at this order. The 
processes which contribute at 0{e^P) are those in which the photon interacts with 
the nucleon's pion cloud. Recent experimental values for the proton polarizabilities 
are @: 

,-(-0.20 ^ in-4f„3 



ap -I- /3p = 13.23 ± 0.86+^;^^ x 10"* fm'^ 



ap -/3p = 10.11 ± 1.74tJ:^2 X 10^4 fm^ , (134) 



where the first error is a combined statistical and systematic error, and the second 
set of errors comes from the dispersion relation analysis used to extract a-p and 
/3p from the experimental data. These values are in good agreement with the xPT 
predictions. 

Meanwhile, the neutron polarizabilities and /3n are difhcult to measure, due 
to the absence of free neutron targets, thus this prediction is not well tested. Hence 
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Fig. 18. The deuteron charge and quadrupole form factors to order eP^ in chiral per- 
turbation theory, with the effects of nucleon structure included in the calculation. The 
experimental data is taken from the compilation of Ref. . Gq is in units of fm^ . 



there has been a considerable effort to extract an and /3n by performing Compton 
scattering on nuclear targets. Coherent Compton scattering on a deuteron target 
has been measured at E-^ = 49 and 69 MeV by the Illinois group |82j and — 
84.2 - 104.5 MeV at Saskatoon g]. 

However, the amplitude for Compton scattering on the deuteron involves mech- 
anisms other than Compton scattering on the individual nucleons inside the deu- 
terium nucleus. Hence, if we want to extract neutron polarizabilities we will need 
a theoretical calculation of Compton scattering on the deuteron in which we can 
make sure that the error in our extraction of an and /?„ from subtracting out these 
"two-body" mechanisms is under control. xPT is an ideal framework for this, since 
the two-body mechanisms that contribute to 7d — > 7d are calculable order-by-order 
in the chiral expansion, and so we can perform a consistent extraction of and 
/3n- The theoretical error is defined as the size of the mistake made by neglecting 
terms of the next order in the expansion. 
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(d) 

Fig. 19. Graphs which contribute to Compton scattering on the deuteron at 0{e^) (a) 



and 0{e^P) (b-d). The sliced and diced blobs are from C'^^^ (c) and (d). Crossed 

graphs are not shown. 



The Compton amplitude we wish to evaluate is (in the 7c? center-of-mass frame) : 



i'M- [ P + ) T^N,,,(k',k) V-Af (P) 



(27r)3 
(27r)6 



2 



^M'(p') r^2Nj^^,^(k',k) v^m(p) 



(135) 



where M (M') is the initial (final) deuteron spin state, and A (A') is the initial (final) 
photon polarization state, and k (k') the initial (final) photon three-momentum, 
which are constrained to |k| = |k'| = w. The amplitude T-yN represents the graphs 
of Fig. 19 and Fig. po| b where the photon interacts with only one nucleon. These 
are the "one-body" mechanisms we actually wish to investigate. Meanwhile, the 
amplitude T^nn represents the graphs of Fig. |20| a where there is an exchanged 
pion between the two nucleons. These are the "two-body" mechanisms we wish to 
calculate in a controlled way and then subtract them off. 

The LO contribution to Compton scattering on the deuteron is shown in Fig. p^(a). 
This graph involves a vertex from c!"^\ and so is O(e^). This contribution is simply 
the Thomson term for scattering on the proton. There is thus no sensitivity to 
either two-body contributions or nucleon polarizabilities at this order. At 0{e^P) 
there are several more graphs with a spectator nucleon (Figs. |9|(b),(c),(d)), as 
well as graphs involving an exchanged pion with leading order vertices (Fig. [2C|(a)) 
and one- loop graphs with a spectator nucleon (Fig. |2^(b)) |84 . Graphs such as 
Fig. p^(b) contain the physics of the proton and neutron polarizabilities at 0{e^P) 
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in xPT. 

The results shown below were generated with the NLO chiral wave function of 
Ref. ||7l|. Fig. ^ shows the results at E-y = 49, 69, and 95 MeV. For comparison 
we have included the calculation at O(e^) in the kernel, where the second contri- 
bution in Eq. (135) is zero, and the single-scattering contribution is given solely by 
Fig. ^(a). At 0{e^P) all contributions to the kernel are fixed in terms of known 
pion and nucleon parameters, so to this order xPT makes predictions for dcuteron 
Compton scattering. 

We also show the O(e^P^) result from work just being completed This cal- 
culation uses the NLO wave function of Ref. and the O(e^P^) single- nucleon 
Compton amplitudes derived in Ref. |85| . Also included are a number of 1 /M cor- 
rections to the two-body piece of the 7NN 7NN amplitude. It turns out that such 
mechanisms are the only corrections to the two-body piece of this amplitude at this 
amplitude. In other words, no two-body counterterms contribute to 7d — > 7d at 
0(e2p2) On the other hand, at ©(e^P^) 

single-nucleon counterterms which shift 
the polarizabilities enter the calculation. But that is an advantage! It is exactly 
these counterterms — or more specifically the neutron counterterms — that we wish 
to extract from the 7d data, since these counterterms serve to fix the value of an, 
/3nj ctp and /3p. In this sense Compton scattering on deuterium at O(e^P^) is anal- 
ogous to the reaction 7d Tr'^d |6j]: an O(e^P^) calculation allows us to test the 
single-nucleon physics which is used to predict the results of coherent scattering 
on deuterium, since there are no undetermined parameters in the two-body mecha- 
nisms that enter to this order in the chiral expansion. For the calculations displayed 
here we have set the nucleon polarizabilities to the values: 



kp = an = 12.0 X 10 
/3p = /3„ = 2.0 X 10 



fm^, 
-4 fm3. 



(136) 
(137) 



In fact the results for this observable are quite sensitive to a-^ and /?n, so ultimately 
one would like to use experimental values for the proton polarizabilities, such as 
those quoted in Eq. (|l34| ), as well as performing a fit to the photon- deitteron 




(b) 



Fig. 20. Graphs which contribute to Compton scattering on the deuteron at 0{e^P). 

Crossed graphs are not shown. 
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Fig. 21. Results of the 0{e'^) (dotted line), 0{e^P) (dashed line), and 0{e^P'^) (solid 
line) calculations for — 49 MeV, 69 MeV and 95 MeV respectively from left to right. 



experimental data in order to constrain the neutron a and (3. This should become 
a more realistic option as data on 7d — > 7d from Lund becomes available . 
The data shown in Fig. includes six points from the Illinois experiment 



at 49 and 69 MeV fSg) and the Saskatoon data at 95 MeV ||3|. Statistical and 
systematic errors have been added in quadrature. It is quite remarkable how well 
the O(e^) calculation reproduces the 49 MeV data. However, since there are sizeable 
higher-order corrections this agreement might be regarded as somewhat fortuitous 
^) . Meanwhile, the agreement of the O(e^P) calculation with the 69 MeV data is 
very good, although only limited conclusions can be drawn. Finally, the results at 
95 MeV disagree markedly with the Saskatoon data, especially at backward angles. 
While this could be indicative that the neutron polarizabilities are starkly incorrect 
tests suggest that modifying a„ and /?„ in order to reproduce the backward-angle 
7d data from Saskatoon would lead to disagreement between our calculations and 
the 69 MeV data. 

In fact, similar conclusions have been reached in potential-model calculations 
p8| , pgj where the results for the differential cross sections are comparable to those 
obtained here using xPT. They are also not dissimilar to those obtained in NN 
EFTs without explicit pions (see Ref. [^ ). 

Many other processes have been computed using Weinberg power-counting. The 
application to 7rd scattering is intriguing, and holds promise of allowing for an 
extraction of the isoscalar ttN scattering length [^3[ |6^. The prediction of the 
threshold amplitude for 7d —> ir'^d was a significant triumph for xPT applied to 
deuterium |Q. Meanwhile, Park, Min, Kubodera, Rho, and collaborators have cal- 
culated many processes, paying particular attention to astrophysical reactions. See 
Ref. g g for summaries of this work, which represents an important application 
of xPT in particular, and effective theories of nuclear physics in general. 



*) In fact, the "Weinberg" power counting used here breaks down at these lower photon energies, 
since it is designed for ui ~ rrin , and does not recover the deuteron Thomson ampUtude as ^ 0. 
Correcting the power counting to remedy this difficulty appears to improve the description of the 
49 MeV data, without significantly modifying the higher-energy results [M. 
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